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ABSTRACT 


Recent  military  applications  have  demanded  photodetectors  with  high  sensitivity, 
high  selectivity  and  multispectral  capability  for  detection  and  identification  of  the  target. 
These  characteristics  have  been  found  in  quantum  well  infrared  photodetectors  (QWIP). 
Driven  by  these  applications,  a  QWIP  photodetector  capable  of  detecting  simultaneously 
infrared  emissions  within  near  infrared  (NIR),  mid  wavelength  infrared  (MWIR)  and  long 
wavelength  infrared  (LWIR)  was  studied,  modeled,  designed  and  characterized.  Using 
the  envelope  function  approximation,  the  mathematical  model  of  the  quantum  phenom¬ 
ena  in  semiconductor  heterostructures  was  derived.  A  computational  tool  was  developed 
to  solve  self-consistently  the  Schodinger-Poisson  equation  using  the  shooting  method,  al¬ 
lowing  the  theoretical  evaluation  of  the  absorption  coefficient.  A  three-color  (NIR, 

MWIR  and  LWIR)  GaAs-based  QWIP  sample  and  a  two-color  (NIR  and  MWIR)  InP- 
based  QWIP  sample  were  designed,  both  comprised  of  stacks  of  uncoupled  wells  for  each 
band  detection.  The  67  layers  of  the  GaAs  sample  was  grown  using  molecular  beam  epi¬ 
taxy  (MBE).  Intersubband  absorption  in  the  sample  was  measured  for  the  MWIR  and 
LWIR  using  Fourier  transform  spectroscopy  (FTIR)  and  the  measured  peak  positions, 
found  at  5.3  pm,  8.7  pm,  and  13.8  pm  are  within  0.3  pm  of  the  theoretical  values,  indi¬ 
cating  that  the  model  accurately  predicts  the  absorption  wavelengths.  A  two-dimensional 
ordered  grating  pattern  was  selected  and  optimized  separately  for  both  MWIR  and  LWIR 
desired  peaks.  Finally  the  photodetector  device  configurations  were  designed  to  permit  to 
the  measurement  of  the  NIR  band  through  photocurrent  spectroscopy  and  performance 
analysis.  The  fabrication  and  characterization  of  the  prototypes  are  a  matter  for  future 
work. 
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EXECUTIVE  SUMMARY 


Recent  military  applications  have  demanded  photodetectors  with  high  sensitivity, 
high  selectivity  and  multispectral  capability.  Both  mecrurium  cadmiun  telluride 
(HgCdTe)  photodiodes  and  quantum  well  infrared  photodetectors  (QWIP)  offer  multi¬ 
color  capability  in  the  medium-wave  IR  (MWIR)  and  long-wave  IR  (LWIR)  spectral  re¬ 
gions.  Each  of  these  technologies  has  its  advantages  and  disadvantages.  However,  the 
possibility  to  achieve  multicolor  detection,  going  from  near  infrared  (NIR)  to  LWIR,  is 
more  likely  to  be  obtained  using  QWIPs. 

Driven  by  those  applications,  and  in  continuity  of  the  work  being  conducted  at  the 
Sensor  Research  Laboratory  (SRL)  at  the  Naval  Postgraduate  School  (NPS),  QWIPs  ca¬ 
pable  to  detect  simultaneously  3  different  IR  bands  within  the  wavelengths  intervals  of 
0.9  -  1.4  pm,  3.8  -  5.0  pm  and  8.0  -  12.0  pm  were  studied,  modeled,  and  designed.  The 
work  was  performed  bounded  by  technological,  financial  and  temporal  limitations,  aim¬ 
ing  the  fabrication  of  proof-of-concept  QWIP  prototypes  to  be  investigated  and  optimized 
in  future  work. 

When  one  semiconductor  layer  (well)  is  sandwiched  between  two  layers  of  larger 
bandgap  material  (barriers),  it  forms  a  quantum  well.  The  potential  profile  is  defined  by 
the  bands  offset,  allowing  the  existence  of  the  quantized  energy  levels,  confined  inside 
the  wells,  while  continuum  states  are  possible  outside.  When  a  quantum  well  photodetec¬ 
tor  is  exposed  to  an  incident  photon  flux,  transitions  between  quantized  energy  levels  can 
occur.  The  photon  flux  can  be  sensed  when  bias  is  applied  to  the  structure  and  the  transi¬ 
tions  allow  the  extraction  of  the  electrons  from  the  well,  to  the  barrier  regions  leading  to  a 
photocurrent. 

The  semiconductor  material  properties  along  with  their  availability  were  consid¬ 
ered  and  two  main  material  groups  were  selected  to  be  used.  GaAs/AlGaAs/InGaAs  sys¬ 
tems  proved  easy  in  growth  and  fabrication  while  InP/AlGaAs/AlInAs  systems  were  ca¬ 
pable  of  meeting  the  requirements  for  detecting  the  desired  NIR  band.  A  configuration  of 


xv 


three  uncoupled  quantum  wells  stacks,  one  stack  for  each  band,  placed  between  highly 
doped  contact  layers,  proved  to  be  more  suitable  to  allow  separate  readouts  of  the  signals. 

To  be  able  to  estimate  the  QWIPs’  absorption  the  quantized  energy  levels  inside 
the  wells  as  well  as  their  respective  wavefunctions  were  computed  self-consistently  solv¬ 
ing  the  Schrodinger-Poisson  equation  numerically  for  the  structures.  The  shooting 
method  was  used  due  to  its  ability  to  handle  any  potential  profile  making  the  design  more 
flexible.  Finally,  a  GaAs-based  sample  for  three-color  detection,  and  a  InP-based  sample 
for  two-color  detection  were  designed,  taking  into  account  the  limitations  on  each  con¬ 
figuration. 

The  GaAs  sample  was  fabricated  using  MBE  and  the  absorption  measurements 
are  presented  in  the  fourth  chapter.  An  FTIR  spectrometer  was  used  to  measure  the  room 
temperature  absorptance  of  the  MWIR  and  LWIR  bands,  and  the  sample  showed  absorp¬ 
tion  in  good  agreement  with  the  theoretical  predictions.  Also  the  heavily  doped  contact 
layers  introduced  a  second  LWIR  absorption  band.  The  measured  peaks,  5.3  pm ,  8.7  pm 
and  13.8  pm,  varied  less  than  0.3  pm  from  the  predicted  values,  mostly  due  to  the  uncer¬ 
tainties  in  the  material  parameters. 

The  optical  coupling  structure  was  than  designed  to  allow  detection  of  normal  in¬ 
cident  radiation.  A  bi-dimensional  ordered  grating  pattern  was  selected  and  optimized 
separately  for  both  MWIR  and  LWIR  desired  peaks.  Constraints  in  mask  fabrication  and 
wafer  processing  forced  the  degradation  of  the  grating  performance,  to  make  it  feasible. 
Finally  the  photodetector  device  configurations  were  designed,  to  permit  to  performance 
analysis.  The  fabrication  and  characterization  of  the  prototypes  are  matter  for  future 
work. 

Despite  the  financial,  temporal  and  technological  limitations,  normally  present  in 
any  research,  the  objective  of  this  thesis  was  successfully  accomplished.  However,  fur¬ 
ther  work  needs  to  be  done  to  finalize  the  devices’  fabrication  and  testing  and  to  improve 
the  mathematical  models  and  the  numerical  algorithms  to  be  able  develop  optimal  detec¬ 
tors. 
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I. 


INTRODUCTION 


Infrared  guided  weapons  have  proved  to  be  indispensable  in  any  modem  battle¬ 
field.  Statistics  and  analysis  have  shown  that  the  capability  to  detect  and  identify  the  en¬ 
emy  platforms’  IR  signatures  as  well  as  the  ability  to  negate  ones  own  are  a  significant 
factors  of  success  in  most  tactical  combat  environments  [1],  In  this  scenario,  the  IR  coun¬ 
termeasures  play  an  important  role.  By  emulating  the  sought  platform  IR  signature  and 
seducing  the  IR  guided  weapon,  it  is  possible  to  reduce  significantly  the  target  risk.  Mod¬ 
em  countermeasures  have  the  capability  to  produce  signatures  very  close  to  the  target’s 
emissions.  Consequently,  to  be  able  to  distinguish  between  countermeasure  and  target, 
the  weapon  seeker  must  have  not  only  high  sensitivity  but  also  high  selectivity  combined 
with  multispectral  detection  capability. 

Systems  that  gather  data  in  separate  IR  spectral  bands  can  discriminate  both  abso¬ 
lute  temperature  and  unique  signatures  of  the  target,  necessary  to  distinguish  it  from  the 
countermeasure.  In  addition,  the  imaging  capability  of  the  focal  plane  arrays  (FPA)  has 
been  explored  in  the  newest  generation  of  seekers,  making  the  weapons  even  harder  to 
counter. 

This  combination  has  been  achieved  until  recently  using  conventional  infrared 
photodetectors  and  extensive  postprocessing  computational  efforts  [2],  Furthermore,  mul¬ 
tispectral  systems  rely  on  complicated  optical  techniques  that  either  disperse  the  optical 
signal  across  multiple  IR  FPAs  or  use  filter  wheels  to  spectrally  discriminate  the  image 
focused  on  a  single  FPA  [3],  In  addition,  beam  splitters,  lenses,  optical  bandpass  filters  as 
well  as  a  complex  alignment  configuration  are  required,  making  these  approaches  expen¬ 
sive  in  size,  complexity  and  cooling  requirements. 

Currently,  both  mecrurium  cadmiun  telluride  (HgCdTe)  photodiodes  and  quantum 
well  infrared  photodetectors  (QWIP)  offer  multicolor  capability  in  the  medium-wave  IR 
(MWIR)  and  long-wave  IR  (LWIR)  spectral  regions.  Each  of  these  technologies  has  its 
advantages  and  disadvantages,  detailed  in  Ref.  [4],  However,  the  possibility  to  achieve 
multicolor  detection,  going  from  near  infrared  (NIR)  to  LWIR,  is  more  likely  to  be  ob¬ 
tained  using  QWIPs. 
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A.  QUANTUM  WELL  INFRARED  PHOTODETECTORS 

The  ability  to  sense  optically  driven  transitions  between  two  quantized  energy 
levels  is  the  origin  of  the  photodetection  in  quantum  well  photodetectors.  Conventional 
interband  optical  absorption  involves  photoexiting  carriers  across  the  bandgap.  Since  the 
incident  photon  energy  is  greater  than  the  bandgap,  both  electrons  and  holes  can  be  cre¬ 
ated  in  the  semiconductor.  In  a  detector,  these  carriers  are  collected,  thereby  producing  a 
photocurrent.  By  controlling  the  bandgap  energy,  selecting  the  desired  material  and  con¬ 
trolling  the  alloy  composition,  the  spectrum  of  the  absorption  can  be  tailored.  [5] 

In  1969  L.  Esaki  and  R.  Tsu  [6]  started  the  investigation  of  the  idea  that  a  periodic 
variation  on  of  the  alloy  composition  or  of  impurity  density  introduced  during  epitaxial 
growth  could  create  a  periodic  potential  or  superlattice,  in  monocrystalline  semiconduc¬ 
tors.  This  structure  is  responsible  for  the  formation  of  subbands  in  both  the  valence  and 
conduction  bands.  In  consequence,  intersubband  optical  absorptions  are  possible  in  quan¬ 
tum  well  photodetectors,  permitting  the  utilization  of  large  bandgap  materials  to  detect 
long  wavelength  infrared  radiation.  Beyond  that,  in  quantum  wells,  both  interband  and  in¬ 
tersubband  absorption  are  possible  to  be  detected  allowing:  monolithically  integrated 
multi  spectral  detection  from  NIR  to  LWIR,  very  high  speed,  and  narrow  and  tunable 
band  detection.  Most  of  those  characteristics,  necessary  for  the  application  in  discussion, 
are  difficult  to  be  achieved  with  conventional  photodetectors  [4], 

Several  studies  and  experiments  were  conducted  since  Esaki’s  paper  confirming 
the  feasibility  of  practical  devices  using  this  concept.  The  first  quantum  well  laser  [7]  was 
reported  in  1975;  the  first  successful  fabrication  of  a  high  electron  mobility  transistor 
(HEMT)  [8]  was  reported  in  1980.  Finally  in  1987  B.  F.  Levine  et  al.  published  the  ex¬ 
perimental  results  of  the  first  working  QWIP  [9],  The  device  was  a  periodic  repetition  of 
AlGaAs/GaAs  and  based  in  conduction  band  bound-to-bound  transitions.  In  the  follow¬ 
ing  years,  various  configurations  have  been  tested  with  great  improvement  in  the  detec¬ 
tion  characteristics,  reaching  the  state  of  art  with  the  first  640x512  pixel  four-band  QWIP 
FPA  camera,  described  in  detail  in  [10], 

In  2002,  a  series  of  investigations  in  dual  color  and  band-tunable  QWIPs  started 
to  be  conducted  at  the  Sensor  Research  Laboratory  (SRL)  at  the  Naval  Postgraduate 
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School  (NPS).  The  results,  published  in  [11],  [12],  [13],  [14],  and  [15],  showed  the  poten¬ 
tial  of  these  devices  in  military  applications.  The  unique  characteristic  of  the  SRL  work  is 
the  inclusion  of  the  NIR  band  in  its  QWIPS,  using  interband  transitions  in  asymmetric 
quantum  wells.  Currently  the  research  has  reached  a  stage  that  demands  more  elaborated 
numerical  models,  devices  capable  of  detecting  normal  incident  IR  radiation  and  more 
than  two  bands  of  military  interest  monolithically  detected. 


B.  PURPOSE  AND  ORGANIZATION  OF  THIS  THESIS 

Driven  by  those  applications,  and  in  continuity  of  the  work  being  conducted  at  the 
SRL,  the  purpose  of  this  thesis  was  to  study,  to  model,  and  to  design  a  QWIP  capable  to 
detect  3  different  IR  bands  within  the  wavelengths  intervals  of  0.9  -  1.4  pm ,  3.8  -  5.0 
pm  and  8.0  - 12.0  pm .  The  work  was  performed  bounded  by  technological,  financial 
and  temporal  limitations,  aiming  at  the  fabrication  of  proof-of-concept  QWIP  prototypes 
to  be  investigated  and  optimized  in  future  work. 

In  order  to  document  the  conducted  research  work,  this  thesis  is  organized  as  fol¬ 
lows.  This  chapter  elucidates  the  motivation  of  this  work  as  well  as  gives  a  brief  introduc¬ 
tion  of  the  research  realm  of  QWIPs.  A  short  description  of  the  organization  of  the  docu¬ 
ment  is  also  presented. 

The  second  chapter  explores  the  physics  of  the  semiconductor  heterostructures 
and  quantum  wells.  The  envelope  function  approximation  is  described  along  with  the 
one-dimensional  potential  solutions  of  the  Schrodinger’s  equation.  The  density  of  states 
and  the  charge  distribution  potential  are  analyzed  providing  a  theoretical  background  to 
understand  the  transitions  between  quantized  energy  levels  within  the  hetereostructures. 
The  possible  quantum  well  transitions  are  analyzed  and  used  to  derive  the  absorption  co¬ 
efficient  in  all  cases  of  interest  of  this  work.  The  chapter  ends  by  providing  the  necessary 
background  information  to  be  employed  in  designing  issues.  The  mathematical  models  of 
the  main  figures  of  merit,  such  as  absorption  spectra,  quantum  efficiency,  responsivity, 
dark  current,  and  detectivity  are  also  discussed. 

The  third  chapter  is  the  core  of  the  document  and  starts  with  the  initial  design 


considerations.  The  basic  device  configuration  is  defined  based  on  the  material  analysis, 
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limitations,  and  constraints,  as  well  as  the  desired  detection  conditions,  light  coupling  and 
restrictions  on  the  bands.  A  specific  numerical  model  that  mathematically  describes  the 
quantum  well  structures,  predicts  the  absorption  line  shapes  and  helps  to  design  band- 
specific  photodetectors  is  developed,  based  on  the  self-consistent  Schrodinger-Poisson 
solutions. 

The  fourth  chapter  presents  the  experimental  results.  Absorptance  measurements 
using  Fourier  transform  infrared  spectroscopy  of  the  molecular  beam  epitaxially  (MBE) 
grown  wafer  are  presented  and  compared  with  the  estimated  values.  The  analysis  and 
conclusions  are  shown  and  the  possible  future  improvement  and  optimization  are  dis¬ 
cussed.  Sequentially,  the  optical  coupling  is  analyzed  and  a  two-dimensional  grating 
structure  is  designed  to  couple  the  desired  IR  bands.  The  final  device  structure,  as  well  as 
several  photodetector  configurations,  are  designed  to  meet  the  standards  and  rules  of  fab¬ 
rication.  The  layout  of  the  masks  is  finalized  in  accordance  with  the  processes  and  budget 
restrictions. 

The  conclusive  chapter  summarizes  the  research  work,  discussing  the  main  as¬ 
sumptions,  along  with  the  technological,  temporal  and  financial  limitations  faced  during 
the  thesis  work  period.  Also  the  lessons  learned  are  addressed.  The  chapter  ends  pointing 
the  scientific  contribution  of  the  work  and  identifying  the  most  important  needs  of  im¬ 
provement  and  continuity  of  the  research. 
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II.  QUANTUM  WELL  INFRARED  PHOTODETECTORS 

PHYSICS 


The  understanding  of  the  physical  concepts  involved  in  IR  photodetection  by 
quantum  well  devices  lies  in  the  quantum  physics  realm,  as  the  name  suggests.  The  abil¬ 
ity  to  control  such  phenomena  requires  the  capability  to  mathematically  model  quantum 
effects  caused  by  semiconductor  structures  and  its  interacting  fields.  This  chapter  pro¬ 
vides  the  theoretical  background  necessary  to  physically  understand  and  practically  de¬ 
sign  QWIPs  through  the  manipulation  of  the  mathematical  models. 

Initially,  the  one-dimensional  potential  present  in  heterostrucures  is  described  by 
Schordinger’s  equation,  along  with  the  concepts  necessary  to  understand  how  quantum 
well  devices  absorb  IR  and  how  to  evaluate  and  predict  such  absorption.  After  that, 
QWIPs  are  described  and  their  figure  of  merit  are  derived  to  permit  predicting  and  ana¬ 
lyzing  their  performances. 


A.  SEMICONDUCTOR  HETEROSTRUCTURES 

One  of  the  most  important  characteristics  of  semiconductors  is  the  existence  of 
forbidden  energy  bands,  responsible  for  separating  the  electrons  in  the  highly  populated 
valence  band  from  those  in  the  poorly  populated  conduction  band.  The  energy  gap  and 
bands  are  determined  by  the  periodic  potential  of  the  crystalline  material.  When  one 
semiconductor  layer  is  grown  on  top  to  the  other,  a  heterojunction  is  formed.  This  can  be 
obtained  using  selected  semiconductors  with  compatible  crystal  structures  and  lattice 
spacings.  Using  modem  techniques  such  as  Molecular  Beam  Epitaxy  (MBE)  and  Metal 
Organic  Chemical  Vapor  Deposition  (MOCVD),  under  carefully  controlled  conditions, 
the  compositional  transition  between  two  materials  can  be  made  almost  perfectly  abrupt 
with  the  heterointerfaces  being  defined  on  a  monolayer  scale.  [16] 

Away  from  the  heterojunction,  the  electrons  are  subjected  to  the  bulk  properties 
of  each  material.  At  the  heterojunction,  the  potential  changes  abruptly  from  one  material 
to  the  other.  Moreover,  a  transfer  of  electrical  charge  occurs  over  the  scale  of  a  few 
atomic  layers  near  the  interface.  This  redistribution  of  charges  creates  an  interface  dipole 
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responsible  for  a  sharp  jump  in  the  electrostatic  potential  that  is  superimposed  to  the  crys¬ 
talline  potential.  On  each  side,  the  electrons  take  on  the  characteristics  of  the  constituent 
semiconductors.  The  relative  positions  of  the  bands,  however,  are  considered,  from  a 
macroscopic  viewpoint,  as  a  discontinuity  or  “offset”  in  the  valence  and  conduction 
bands,  that  depends  uniquely  on  the  pair  of  materials  forming  the  heterojunction.  [5] 

In  terms  of  the  band  discontinuities  three  types  of  structures  are  possible,  classi¬ 
fied  into  four  groups  [16]:  type  I,  type  II-misaligned,  type  II-staggered,  and  type  III. 

Type  I  occurs  when  the  entire  gap  of  the  small  bandgap  semiconductor  resides  within  the 
conduction  and  valence  bands  of  the  large  bandgap  material.  Type  II,  alternatively,  oc¬ 
curs  when  one  (staggered)  or  both  (misaligned)  band  offsets  are  larger  than  the  difference 
between  the  semiconductors  bandgaps.  Finally,  type  III  occurs  when  the  heterojunction  is 
formed  by  one  semiconductor  with  positive  bandgap  and  the  other  with  a  negative  band- 
gap.  The  last  behaves  like  semimetals,  since  the  conduction  and  valence  bands  overlap. 

Taking  into  account  the  bulk  properties  of  the  semiconductors,  it  is  possible  to 
create  virtually  any  potential  profile  along  the  crystal  growth  direction,  by  depositing  se¬ 
quences  of  different  materials.  Limitations  such  as  the  material  availability  and  their  ca¬ 
pability  to  be  grown  epitaxially  as  well  as  the  difference  in  lattice  constants  must  be  con¬ 
sidered  in  any  design.  In  the  rest  of  the  chapter  a  detailed  discussion  on  the  electron  states 
in  heterostructures  is  presented. 

1.  Envelope  Function  Approximation 

The  electron  wavefunctions  in  a  heterostructure  can  be  described  using  the  enve¬ 
lope  function  approximation,  as  described  in  Ref.  [5], 

A  semiconductor  heterostructure  can  be  considered  as  possessing  a  periodic  po¬ 
tential  Vc(r),  due  to  the  crystal  lattice,  and  a  slowly  varying  electrostatic  potential  Ve(r) 
due  to  different  layers  of  the  heterojunction.  Here  r  represents  the  position  vector.  The 
Hamiltonian  for  an  electron,  given  this  combination  of  potentials  can  be  written  as  [5]: 

H  = -^— +  Vc(r)  +  Ve(r) .  (2.1) 

2m0 
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The  momentum  operator  p  is  described  as 


d_  d_  d_ 

dx  ’  dy  ’  dz 


(2.2) 


where  h  is  the  Plank  constant  h  divided  by  2 n  and  i  is  the  imaginary  unit.  Along  with 
some  very  demanding  algebra,  some  physically  consistent  approximations  [5]  can  be 
made  under  the  assumption  that  the  electron  in  the  nth  band  is  therefore  a  particle  pos¬ 
sessing  effective  mass  associated  with  the  nth  band  and  subjected  to  a  potential  Ve(r). 
Furthermore,  considering  that  influence  of  the  periodic  potential  Vc(r)  is  taken  into  ac¬ 
count  through  the  effective  mass  m  *n  and  the  energy  En0  of  the  nth  band  at  k  =  0  (k  be¬ 
ing  the  crystal  momentum),  the  Schrodinger’s  equation  for  the  envelope  function  can  be 
written  as  [5]: 

H'¥n(r)  =  (E-En0)'¥n(r)-  (2-3) 


Here,  His  the  Hamiltonian  described  by  Equation  (2.1)  (replacing  mo  by  the  effective 
mass)  and  E  represents  the  energy  eingenstates  of  the  wavefunctionvPn .  The  full  wave- 
function  can  be  written  as  [5] 

T„=^(rK0(r),  (2.4) 

where  y/n  is  the  envelope  function,  and  un0  is  the  Bloch  function  possessing  the  periodic¬ 
ity  of  the  lattice.  When  the  electrostatic  potential  is  zero,  the  solution  is  [5] 

'f'"S^2C*r""0<r)’  (2'5) 


where  £2  is  the  crystal  volume. 

Further  simplifications  are  possible,  taking  into  account  the  one-dimensional 
characteristic  of  the  heterostructure  potential  Ve(r)  and  the  effective  mass  along  the 
growth  axis  z.  Consequently  the  envelope  function  can  be  written  in  the  form  [5] 

¥i  r)  =  ^(z)e'Kp.  (2.6) 
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It  is  important  to  distinguish  between  r  (that  is  the  3-D  position  vector)  and  p  (that  is  the 
2-D  position  vector  in  the  x-y  plane).  Also,  while  k  represents  the  3-D  crystal  momentum, 
K  is  the  2-D  representation  of  the  momentum  in  the  x-y  plane. 

The  z  component  of  the  envelope  function  satisfies  [17] 

~  \(  \  7  2  Y(z)  +  V(z)yf(z)  =  Ey/(z).  (2.7) 

2m  *  (z)  dz 

To  obtain  the  stationary  states  of  an  electron  in  the  heterostructure,  it  is  necessary 
to  solve  Equation  (2.7)  in  each  layer,  considering  the  properties  of  each  component  sepa¬ 
rately  and  to  join  the  solutions  at  the  interfaces  [17].  First,  the  wavefunctions  must  be 
continuous  at  the  interface.  Mathematically, 

W(zi)  =  W(zi+)>  (2-8) 

where  z,  represents  the  z  coordinate  of  the  /th  interface  and  the  (±)  sign  refers  to  the  mate¬ 
rials  on  the  right  and  left  of  the  heterojunction.  Second,  the  solution  must  conserve  the 
probability  flux  which  crosses  the  interface.  This  requirement  is  satisfied  when  [18] 

— fr(V)  =  -rf^/+)>  (2-9) 

m  dz  m  dz 

Finally,  the  probability  interpretation  of  the  wave  function  [17]  requires  the  nor¬ 
malization  condition  for  the  bound  states  potential.  Mathematically, 

J  y/  *  (z)^(z)  dz  - 1 ,  (2.10) 

where  (*)  represents  the  complex  conjugate. 

Some  important  constrains  of  this  formalism  must  be  mentioned  [18].  First,  the 
effective  mass  was  assumed  to  not  depend  on  k,  meaning  a  perfect  parabolic  band.  When 
the  discontinuities  in  V(z)  are  significant,  a  correction  accounting  for  band  nonparabolic- 
ity  must  be  applied.  This  can  be  done  using  the  expression  [19] 

>»♦(*)  =  >»,  *(z)(l  +  y B(E  -  F,(z>)),  (2.11) 
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where  mp  *  (z)  is  the  “parabolic”  band  effective  mass  profile  and  [3  is  given  by  [19]: 


P  = 


1- 


m*(z ) 
m0 


y 

J 


1 

w 


(2.12) 


Here,  mo  is  the  electron  rest  mass  and  Eg(z)  is  the  bandgap  profile  along  the  heterostruc¬ 
ture.  Second,  all  parameters  characterizing  the  band  structure  are  temperature  dependent 
This  dependence  can  be  modeled  by  the  Varshini  equation  [19]: 

a  T2 

=  (2.13) 

1  Ptmp 


where  E  0  is  the  bandgap  at  T  -  0 ,  and  atmp  and  ft  are  adjustable  Varshini  parame¬ 
ters  that  can  be  found  for  several  semiconductors  compounds  in  [20]. 

Finally,  the  strain  placed  on  the  material  by  growing  layers  with  different  lattice 
constants  on  top  of  each  other  also  affects  the  band  structure.  This  allows  an  additional 
degree  of  freedom  in  designing  heterostructure  devices.  For  lattice-mismatched  epitaxial 
growth  on  a  (001)  substrate,  the  strain  tensor  is  given  by  [21] 


a,  1 

£  =  £  =  — -1 

xx  yy 


£  =  -- 


2  a 


12 


(2.14) 


where  as  and  aL  stand  for  the  lattice  constants  of  the  substrate  and  the  strained  layer,  re¬ 
spectively,  and  cj2  and  cn  are  elastic  constants  for  crystals  with  cubic  cells  under  small 
deformations.  All  off-diagonal  terms  are  zero  in  this  case.  Therefore,  applying  the  defor¬ 
mation  potential  theory,  it  is  possible  to  compute  the  bandgap  shift  [21].  For  III-V  mate¬ 
rials,  due  to  the  nature  of  the  atomic  bonding,  the  bandgap  increases  for  compressive 
strain  and  the  change  is  given  by  [20] 

5Eg  =  (au  +  ac )(exx  +  eyy  +  ezz) ,  (2.15) 


where  av  and  ac  are  the  empirical  deformation  potentials,  also  given  in  [20], 
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Rigorously,  the  assumptions  made  until  now  cannot  be  applied  to  the  valence 
bands  that  are  degenerate  with  respect  to  the  heavy  and  light  hole  bands.  A  formal  treat¬ 
ment  for  the  valence  subbands  can  be  found  in  [21]  and  the  envelope  function  approxima¬ 
tion  using  a  4x4  k-p  Hamiltonian  approximation  can  be  compiled  as 


P  +  Q 

L  M 

id 

P-Q  o 

0  P-Q 

0 

-L' 

0 

M 
-L 
P  +  Q 


^rl(k) 

^2(k) 

^4(k) 


=  £(k) 


^‘(k) 

^2(k) 

^3(k) 

^4(k) 


(2.16) 


where  \i/vi  is  the  envelope  wavefunction  of  the  holes  under  the  spin  symmetry  specified 
as  follows:  v\  and  v4  represents  the  1 3/2, ±3/2)  heavy  hole  spin  states,  and  v2  and 
v3 represents  the  1 3/2, ±1/2)  light  hole  spin  states.  The  symbol  (1 )  stands  for  adjoint, 
and  the  matrix  elements  are  given  by  [21]: 


P  = 


2  m0 


C K2+K 2)- 


h2  d 
2  m0  dz 


Vi(z)-1-  +  Ve(z) > 
dz 


(2.17) 


Q  = 


2  m0 


(K2  +  K2)  + 


tf  d 

2m0  dz 


ri(z)^+Re(z), 

dz 


(2.18) 


fi24 3  .,  , 

L  =  —~ - (kx~lky) 

2  mn 


,  ,  d  d  .  . 

r3(z)— +— r3(z) 

dz  dz 


(2.19) 


and 

m =PP-[(r2(z)-r,(z)](kx -ikyf.  (2.20) 

4w0 

All  parameters  in  the  previous  expressions  were  defined  before  except  for  the  “Luttinger 
parameters”,  yi ,  that  can  be  found,  for  several  III-V  compound  semiconductors  and  their 

alloys,  in  [20].  It  is  important  to  mention  that  the  z  dependence  of  yi  and  V  is  due  the  ma¬ 
terial  change  along  the  growth  axis. 

For  practical  applications,  further  simplifications  are  acceptable.  The  most  com¬ 
mon  is  to  consider  the  valence  band  represented  by  two  parabolic  bands,  one  for  the 
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heavy  holes  with  effective  mass  mhh  *  <  0 ,  and  other  for  the  light  holes  with  effective 
mass  mlh*  <  0 .  The  well  potential  is  then  a  confining  potential  for  holes  and  the  same 
model  used  for  electrons  in  the  conduction  band  becomes  applicable. 

In  a  more  practical  sense,  a  heterostructure  can  be  grown  to  create  potential  wells, 
with  electrooptical  proprieties  that  can  be  customized  for  a  wide  variety  of  applications. 
This  practice  is  commonly  called  “quantum  engineering”  or  “bandgap  engineering”. 


2.  Quantum  Wells 

When  one  thin  semiconductor  layer  (well)  is  sandwiched  between  two  layers  of 
larger  bandgap  material  (barriers),  it  forms  a  quantum  well.  As  mentioned  in  the  previous 
section,  the  potential  profile  is  defined  by  the  bands  offsets,  allowing  the  existence  of  the 
quantized  energy  levels,  confined  inside  the  wells,  while  continuum  states  are  possible 
outside.  Such  structures  allow  the  exploration  of  quantum  effects  that  have  become  very 
useful  in  optoelectronic  devices.  [5] 

A  schematic  representation  of  the  electrostatic  potential  of  a  structure  comprised 
of  multiple  quantum  wells  is  shown  in  Figure  1 .  The  barriers  and  wells  are  highlighted 
with  bandgaps,  EgB  and  EgW,  respectively.  The  other  parameters  in  Figure  1  include  va¬ 
lence  and  conduction  band  offsets,  VBO  and  CBO,  respectively;  and  the  valence  and  con¬ 
duction  bands  confined  energy  levels  Evi  and  Ea,  respectively. 


4  V(z) 


Conduction  Band 


A  A 


CBO 


Ec2 

Ed 


-gB 


zgw 


jvBO 


Valence  Band 


Ec2 

Ed 


"i 


Evi  f 

Ev2 


Evi 

Ev2 


Ec2 

Ed 


Evi 

Ev2 


Continuum 


Ec2 

Ed 


Evi 
Ev2 


Continuum 


Figure  1  Schematic  diagram  of  a  heterostructure  potential  profile  Ve(z )  (dark  continu¬ 
ous  line). 
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The  ability  to  explore  this  structure  as  a  device  requires  knowledge  of  the  in¬ 
volved  physical  phenomena  through  mathematical  models.  The  behavior  of  the  carriers 
inside  the  heterostructure  can  be  described  by  the  Schrodinger’s  equation  along  with  the 
appropriate  boundary  conditions.  Complementarily,  the  concepts  of:  density  of  states; 
subbands  population;  charge  distribution  potential;  transition  rates;  and  effects  of  external 
fields;  among  others,  have  to  be  elaborated  in  order  to  allow  practical  comprehension  and 
required  control  in  designing  procedures. 

Therefore,  the  efforts  from  now  on  are  focused  in  describing  mathematically  the 
most  important  characteristics  of  quantum  wells,  considering  their  suitability  to  act  as 
photodetectors. 


3.  Schrodinger’s  Equation  Solutions 

For  the  sake  of  simplicity,  we  consider,  initially,  an  electron  in  a  conduction  band. 
A  schematic  diagram  of  an  TV-layer  multiple  quantum  well  structure  is  shown  in  Figure  2. 


4E 


Figure  2  Conduction  Band  schematics  diagram  of  an  TV-layer  structure.  The  well  width 


is  Lw . 


The  general  solution  of  the  Shrodinger  equation  for  the  nth  layer  can  be  written  as 


lif  =  A 

t  n  i 


,e'v  +  Be 


ik„z 


(2.21) 


where 


K(E)  = 


/2»c  *(£-*;) 

v  n2 


(2.22) 
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The  coefficients,  An  and  Bn,  can  be  evaluated  applying  the  boundary  conditions  described 
by  Equations  (2.8)  and  (2.9),  leading  to  the  matrix  form 


'A 

=  M 

An+\ 

A. 

n 

_Bn+\  _ 

(2.23) 


where 


M 


n 


2 


(i  +  r„yft-ry 

(l-Vny(KM 


{\-Yn)ei(kMz’ 
(1  +  yn)e-iiK-'-K)z’ 


(2.24) 


and 


Yn 


_mn 


Vk 1 


mnAK 


(2.25) 


Here,  if  the  structure  contains  N  layers,  the  coefficients  of  the  first  and  last  layers  are  re¬ 
lated  by 


A 

A 


(2.26) 


where  the  interior  product  turns  out  a  4-by-4  matrix  M  .  The  wavefunction  normalization 
condition  expressed  in  Equation  (2.10)  further  demands  that  i//(z)  —>  0  and 

d\j/{z)/dz  — >  0 ,  as  z  — >  ±oo  for  bound  states.  Therefore,  since  the  two  outer  layers  are 
barriers,  and  for  E  <  Vn ,  the  coefficients  Bj  and  AN  must  be  zero.  Thus,  Equation  (2.26) 
can  be  rewritten  as 


1 

a 

mn 1 

1 - 

o 

1 _ 

1 

o 

_ 1 

_m2 1 

- 1 

<N 

7 
_ 1 

(2.27) 

To  evaluate  the  confined  energy  states  it  is  necessary  to  solve,  numerically,  the  equation 

m22(E)  =  0.  (2.28) 


Finally,  the  wavefonctions  can  be  obtained  making  BN  =  1 ,  for  example,  and  solving  se¬ 
quentially,  for  each  layer,  Equation  (2.23). 
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For  E  >V0,  the  structure  admits  delocalized  states  where  the  electron  can  take  on 

any  value  of  positive  energy,  forming  a  continuum  state.  Since  the  states  are  not  bound 
any  more,  the  normalization  of  the  wavefunction  becomes  a  problem.  A  simple  technique 
to  overcome  that  is  described  later  in  this  chapter. 

This  approach  is  called  the  Transfer  Matrix  Method  [22],  Despite  its  simplicity,  it 
has  proven  to  be  a  powerful  tool  to  handle  structures  with  multiple  layers.  The  confined 
states  as  well  as  their  respective  wavefunctions  are  depicted  in  Figure  3  for  a  single-well 
structure. 


Figure  3  Confined  energy  states  and  respective  wavefunctions  in  a  single  square  well. 

It  is  proved  [17]  that,  for  symmetric  structures  like  the  well  in  the  Figure  3,  the  wavefunc¬ 
tions  are  mutually  orthogonal,  in  the  sense  that 

| ¥n(z)*¥m(z)  dz  =  0  (2.29) 

whenever  n^m  .  This  result  is  a  key  point  to  establish  the  electron  transition  rules. 

When  the  structure  barrier/well/barrier  is  periodically  repeated,  and  the  carrier 
de  Broglie  wavelength  is  of  the  order  of  the  barrier  thickness,  an  energy  miniband  is 
formed.  Thus,  the  wavefunctions  of  individual  wells  tend  to  overlap  due  to  tunneling. 
Such  configurations  are  interesting  to  broaden  the  absorption  spectra  as  well  as  allowing 

photovoltaic  detection  [23].  Nevertheless,  when  the  barriers  are  thick  (greater  than  200A, 
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in  most  cases),  the  wells  in  the  superlattice  are  uncoupled  and  the  structure  exhibits  the 
individual  wells  with  discrete  energy  levels  [21].  This  is  the  case  of  interest  in  this  work 
and  all  treatment  onward  considers  this  configuration. 


4.  Density  of  States  and  Subband  Populations 

To  be  able  to  map  the  distribution  of  energy  and  momentum  of  electrons  in  a 
quantum  well,  it  is  necessary,  first  of  all,  to  determine  the  density  of  states.  Using  a  two- 
dimensional  electron  gas  model  due  to  the  characteristic  of  the  structure,  the  density  of 
states  of  a  single  subband  per  unit  area  ( J  'm-2 )  is  given  by  [18]: 

p2D(E)  =  - (2.30) 
nti 

If  there  are  n  confined  states  within  the  quantum  well,  then  the  density  of  states  at 
a  particular  energy  is  the  sum  over  all  subbands  below  that  point  and  can  be  written  as 

p2D(E)  =  fj^-@(E-Ei)  (2.31) 

i= 1  ft" 

where  ©  is  the  step  unit  function. 

The  surface  density  of  carriers  in  a  subband  is  given  by  [18] 

n„=\  hh  fFD(E)p2D(E)dE,  (2.32) 

J subband 

where  fFD  (Fermi  Dirac  statistics)  is  the  probability  of  occupation  [19] 


fFD(E) 


1  +  e 


E-Ep 

kRT 


(2.33) 


Here,  £>  represents  a  quasi-Fermi  level  which  describes  the  carrier  population  within  a 
subband.  Evaluation  of  Equation  (2.32)  gives  [5] 


n,„  - 


me  *kBT 
nti2 


In 


1  +  e 


ef-e„ 

k„T 


(2.34) 
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The  total  electron  density  in  a  system  is  obtained  by  summing  the  densities  of  all 
subbbands. 

In  practice,  a  good  approximation  is  to  deduce  the  number  of  carriers  directly 
from  the  doping  densities  in  the  well  layer  or  the  surrounding  barrier  layers.  This  allows 
the  evaluation  of  the  quasi-Fermi  energy,  the  only  unknown  in  Equation  (2.34). 

This  concept  is  quite  important  and  necessary  to  determine  the  band  bending 
caused  by  potential  due  to  the  charge  distribution. 

5.  Charge  Distribution  Potential 

In  practical  devices,  the  doping  concentration  can  be  high  enough  to  give  rise  to  a 
significant  additional  potential  on  top  of  the  usual  band-offsets.  In  this  case,  it  becomes 
necessary  to  solve  the  electrostatics  that  describes  the  system.  In  quantum  well  structures 
it  is  reasonable  to  assume  the  charge  density  given  by  the  doping  density  of  the  well 
and/or  surrounding  layers.  In  a  n-doped  quantum  well  structure,  the  electrons  are  trapped 
in  the  well  layer,  where  the  potential  energy  is  the  lowest,  making  this  region  negatively 
charged.  Consequently  the  left  and  right  barriers  become  positively  charged.  The  addi¬ 
tional  potential,  V  (z) ,  arising  form  this  charge  distribution,  p ,  can  be  derived  using 

Poisson’s  Equation  [24] 

^2Vp=—,  (2.35) 

e 

where  e  is  the  permittivity  of  the  material.  The  solution  can  be  obtained  determining  the 
electric  field  strength.  The  potential  and  the  electric  field  are  related  by 

z 

Vp(z)  =  -jE-dz.  (2.36) 

Here,  the  one-dimensional  characteristic  of  the  bandedge  potential  imposes  the  one¬ 
dimensional  charge  distribution. 

In  a  doped  semiconductor,  there  are  basically  two  contributions  to  the  charge  den¬ 
sity,  the  ionized  impurities  and  the  free  carriers.  While  the  former  can  be  obtained  from 
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the  doping  density  profile,  d  (z) ,  the  latter  can  be  computed  from  the  probability  distribu¬ 
tions  of  the  carriers  at  the  heterostructure,  y/i  *(z)tpi(z) .  Thus,  the  charge  density  can  be 
expressed  as 

p(z)  =  q[Ny/*(z)y/(z)-d(z)\,  (2.37) 


where  q  is  the  charge  on  the  extrinsic  carriers  and  N  is  the  total  number  of  carriers.  Fur¬ 
thermore,  the  volumetric  charge  density  p{z)  can  be  considered  as  an  surface  charge 
density  a(z)  with  infinitesimal  thickness  Sz .  This  can  be  expressed  mathematically  as: 

o{z)8z  =  p(z).  (2.38) 


When  more  than  one  subband  is  populated,  then  the  contribution  to  the  charge  density 
must  be  summed  over  the  relevant  subbands,  resulting  in 


cr(z)  =  q 


YJNli//l*(z)y/i(z)-d(z) 

Z=1 


Sz . 


(2.39) 


Since  in  most  of  the  cases  the  dimensions  in  the  x-y  plane  are  much  greater  than 
the  z  direction,  the  x-y  plane  can  be  considered  an  infinite  uniform  charge  distribution. 
Thus,  the  electric  field  at  z  can  be  written  as  [18] 

E(z)  =  X  ^r^sisn(z  ~  z')  •  (2.40) 

z'=-oo 


In  Equation  (2.40),  the  summation  represents  the  contribution  of  each  surface  density 
slice  to  the  electric  field  at  the  position  z.  The  function  “sign”  is  defined  as 


sign(z) 


J+l  forz>0 
[-1  forz<0, 


(2.41) 


and  was  introduced  to  account  on  charge  neutrality  along  the  structure. 

At  this  point,  all  parameters  are  in  place  to  allow  the  evaluation  of  V  (z) ,  given 

by  Equation  (2.36).  However  it  should  be  noticed  that  the  wavefunctions  are  needed  to 
obtain  the  density  of  electrons.  This  recursive  problem  can  be  solved  using  the  so  called 
“Self-Consistent  Schrodinger-Poisson  Solution.”  This  method  is  detailed  described  in  the 
next  chapter. 
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6.  Transitions 

Transitions  between  different  states  in  a  quantum  well  can  occur  in  many  ways. 
Nonetheless,  two  kinds  of  transitions,  caused  by  a  time-varying  electromagnetic  perturba¬ 
tion,  are  of  interest  of  this  work:  an  interband  transition,  in  which  an  electron  in  the  va¬ 
lence  subband  can  be  exited  to  the  conduction  subband,  and  an  intersubband  transition,  in 
which  an  electron  moves  from  one  subband  to  another  while  remaining  in  the  same  band 
(conduction  or  valence).  In  both  cases,  depending  on  the  structure  characteristic,  transi¬ 
tions  can  occur  from  a  confined  state  to  a  continuum  state.  All  four  possibilities  are  dis¬ 
cussed  in  this  section. 


In  the  general  sense,  the  transition  rate  from  initial  a  state  represented  by  vPiK.  to 
final  states  represented  by  'f  /K  ,  due  to  an  interaction  potential,  Vp,  can  be  calculated  us¬ 
ing  Fermi’s  golden  rule  [5] 


(2.42) 


In  Equation  (2.42),  hco  is  the  incident  photon  energy,  and  the  delta  function  accounts  for 
conservation  of  energy.  The  interaction  potential,  under  the  dipole  approximation,  can  be 
written  as  [25] 


Vp  = 


m„ 


in 


y/2 


v2  e0nrcoc  j 


epe 


(2.43) 


where  q  is  the  electron  charge,  I  is  the  incident  photon  flux,  c  is  the  speed  of  light,  nr  is 
the  index  of  refraction,  co  is  the  angular  frequency  of  the  incident  photon  flux,  e  is  the 
unit  vector  of  the  photon  polarization,  and  p  is  the  vector  form  of  the  momentum  opera¬ 
tor.  In  the  exponential  term  s  represents  the  wavevector  of  Vp.  Since  the  main  interest  is 
the  interaction  with  light  and  the  photon  momentum  is  negligible  compared  to  the  elec¬ 
tron  momentum,  the  exponential  term  can  be  dropped  out. 

The  term  in  the  Dirac  notation  in  Equation  (2.42)  is  called  electric  dipole  matrix 
element.  Combining  Equations  (2.4)  and  (2.6)  to  represent  the  wavefunctions  x¥n ,  the 
matrix  element  can  be  expanded,  resulting  in 
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(2.44) 


)  =  ^(^/^)"/o(>-yX'T>,K,  (Z>", »(>•)«*'  ')  • 

Considering  the  fact  that  the  envelope  functions  (related  to  the  band-offset/electrostatic 
potential)  vary  slowly  compared  to  the  Bloch  functions  (related  to  the  lattice  potential) 

and  the  initial  and  final  states  have  the  same  momentum  (K/  =  K;  j ,  Equation  (2.44)  can 

be  rewritten  as  [25] 

(^/K \v)~(uf„(r)\Vp\uK(r)){Vf(z)\V,(z)) 

+ (r ,  (z)  Y,  I V,  (z))  |  u,„(r)) . 

This  expression  can  be  further  simplified  considering  the  specific  aspects  of  each  kind  of 
transition. 


a.  Interband  Transitions 

In  this  case,  the  transitions  occur  between  the  valence  and  conduction 
bands.  Given  that  K;  and  K  f  are  the  same  and  are  far  from  the  edge  of  the  Brillouin 

zone  (otherwise  the  envelope  function  approximation  is  not  valid),  the  Bloch  functions 
are  orthogonal.  This  cancels  out  the  second  term  in  the  right  hand  side  of  Equation  (2.45). 
The  expansion  of  the  interaction  potential,  in  the  remainder,  results  in 


m„ 


in 


y/2 


K2£0nrcoc  j 


w/o(r)|e-p|w.0(r)}(^/(z)|^.(z)}.  (2.46) 


The  insertion  of  Equation  (2.46)  into  Equation  (2.42)  results  in  the  rate  of  transitions  per 
state,  given  by 


W  = 


nq 


(m*)  £0nrcoc 


Vo 


(r)|e-p|w,.0(r) 


x 


Wj(z)\w,(z) 


f  n2K2 

- - +  Ef-Ei-hco 

K2mr*  f  ' 


(2.47) 


where  Et  and  Ef  are  the  initial  and  the  final  confined  energy  levels. 
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The  extra  term  in  the  argument  of  the  delta  function  accounts  for  the  en¬ 
ergy  displacement  due  to  the  value  of  K  (conduction  and  valence  bands  are  assumed  pa¬ 
rabolas  with  opposite  concavity).  The  reduced  mass  mr*  is  defined  as 


1 

mr* 


1  1 

■  +  ■ 


me  *  mh  * 


(2.48) 


where  me*  and  mi*  are  the  effective  masses  of  the  electron  and  hole,  respectively. 

At  this  point  it  is  necessary  to  sum  over  all  states  (K).  Considering  the  K- 
space  a  continuous  space,  and  taking  into  account  all  occupied  states  in  the  valence  band 
and  the  correspondent  available  states  in  the  conduction  band  through  the  occupation  sta¬ 
tistics,  it  is  possible  to  express  the  transition  rate  as 


W  = 


71  q 


(jn*f  e0nra)c 


f  o' 
f 


(r)|e-p|w,.0(r)}  (wf(z)  yrt(z) 


x 


(2  n) 


-d2  K 


2^2 


+  Ef-  E.  -  ha) 

2m..  *  f  1 


(2.49) 


Here,  the  term  out  of  parenthesis  inside  the  integral  is  the  2-D  joint  density  of  states  in¬ 
volved  in  the  transitions,  and  f\FDv  (  £) )  and  [1  -fFDc(Ef)\  represent  the  Fermi-Dirac 

probability  of  occupancy  of  the  valence  band  and  probability  of  unocupancy  of  the  con¬ 
duction  band,  respectively.  The  integral  in  Equation  (2.49)  is  evaluated  substituting 


E  = 


n2x2 

2  mr 


(2.50) 


considering  the  Fermi-Dirac  statistics  constant  for  the  K  values  involved  and  using  the 
properties  of  the  delta  function.  The  result  is 


W  = 


q 


IA 


(m*)  £0nrcoc 


i /0(r)  |e  •  p|  w;0(r))|2 1  (wf(z)  \  yrt  (z) 


X- 


m . 


'-©(h(D-(Ef-Ei))fiFDv(Ei)[l-ffFDc(Ef)]. 


(2.51) 


The  modulus  square  of  the  Bloch  functions  matrix  element  term  in  Equation  (2.51)  must 
be  expanded  to  obtain  the  strength  of  the  transition.  For  practical  reasons  [5],  only  normal 
incidence  (TE  mode)  is  considered.  Thus 
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sPvcf  =  (M/o(x)|^|w,.oW)c°s(6>)  +  (M/o(j)|^|M!.o(j))sin(^) 


(2.52) 


where  pvc  is  the  Kane  matrix  element  in  the  bulk  material,  S  stands  for  selection  rule  due 

to  the  type  of  transition,  and  6  is  the  angle  between  the  e  and  the  x  axis.  Here,  for  polar¬ 
ized  light  one  can  select,  for  example  e  =  ex  and  drop  the  y  term  out.  In  the  general  case, 

it  is  necessary  to  take  the  average  incident  direction  of  e .  The  electron  momentum  x  and 
y  components  can  be  considered  the  same  on  average,  due  to  the  heterostructure  charac¬ 
teristics.  Hence,  Equation  (2.52)  can  be  rewritten  as 


e-Sp  vc  = 


U 


/O' 


0)/*ko0))  |cos(0)  +  sin(#)|2 


(2.53) 


Moreover,  the  average  in  6  for  the  second  term  in  the  right  hand  side  of  Equation  (2.53) 
is  one.  The  matrix  element  left  can  be  evaluated  in  terms  of  Kane  energy,  Ep  (eV),  and  for 
transitions  close  to  the  subband  edges,  where  K  ~  0 ,  it  can  be  approximated  as 


(m/o(x)/Jm,.0(x)) 


(2.54) 


where  Ep  is  unanimously  accepted  in  the  literature  as  varying  between  17  and  25  eV  for 
most  of  the  III-V  binaries.  An  extensive  collection  of  this  parameter  is  given  in  [20]. 

Finally  the  absorption  coefficient  is  defined  as  the  rate  of  transitions  per 
unit  volume  divided  by  the  incident  photon  flux,  and  can  be  expressed  from  the  above 
derivations,  for  TE  incident  light  with  energy  cofi ,  as 


7W 


( m  *)2  £0nrch(h(O )  L 


±E 
2  p 


l//f(z)\l//i(z) 


(2.55) 


x 


e(ha>-(Ef-E,  )[!-//“(£, )]. 


The  selections  rules  S  are  derived  from  the  wavefunctions  that  satisfy  Equation  (2.16)  and 
can  be  thought  as  a  correction  factor  in  the  valence  band  simplification  made  here.  For 
TE  configuration  [24], 


S,f= iA/« 
=  1  A/2, 


(2.56) 
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where  the  indexes  Ih  and  hh  stand  for  light  and  heavy  holes,  respectively. 

Especially  important  is  to  notice  that  the  overlap  of  the  envelope  function, 
(t//f(z)\i//i(z)Xj ,  in  Equation  (2.55)  imposes  selection  rules  as  well.  Due  to  the  orthogonal¬ 
ity  of  the  envelope  functions  in  symmetric  wells,  transitions  between  subbands  with  dif¬ 
ferent  quantum  numbers  are  forbidden.  In  most  of  the  practical  cases,  transitions  between 
the  first  valence  subband  and  the  second  conduction  subband  are  desired.  In  this  case  the 
symmetry  of  the  quantum  well  must  be  broken  and  it  is  commonly  done  using  step  wells. 

Further  simplification  is  possible  when  undoped  wells  are  used  along  with 
low  temperatures.  In  this  case,  the  probability  of  occupied  states  can  be  considered  to  be 
one  in  the  valence  band  and  zero  in  the  conduction  band.  In  this  sense,  Equation  (2.55) 
reduces  to 


(m*f  £0nrch(hco )  ’  2  f 


(2.57) 


X 


V f(.z)\Vi(.zj)  ©( hco-(Ef-E ;)). 


Afterward,  the  absorption  coefficient  can  be  evaluated,  applying  Equation  (2.57)  for  the 
case  of  heavy  and  light  holes. 

Interband  absorptions  are  experimentally  explored  in  bulk  material  and 
widely  used  in  practical  devices.  In  the  case  of  quantum  wells,  interband  absorptions  are 
interesting  mainly  in  optical  modulators  [26]  or  multiband  detection  [11],  Much  more  ex¬ 
plored  in  QWIPs  are  the  intersubband  transitions  that  allow  long  wavelength  detections 
with  quite  large  bandgap  materials. 


b.  Intersubband  Transitions 

Electronic  transitions  between  subbands  within  the  same  band  occur  in 
doped  quantum  wells.  Valence  band  transitions  require  p-type  doped  material  while  con¬ 
duction  band  transitions  require  n-type  doped  semiconductors.  Due  to  the  simplicity  of 
the  conduction  band  compared  to  the  valence  band,  the  conduction-band  configuration  is 
more  practically  explored.  That  is  the  case  of  interest  in  this  work. 
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Starting  with  Equation  (2.45),  the  first  term  in  the  right  hand  side  van¬ 
ishes,  due  to  the  orthogonality  of  the  envelope  functions.  In  addition,  the  Bloch  overlap, 
in  the  same  term,  goes  to  one  since  the  Bloch  functions  are  normalized  to  unity.  The  ex¬ 
pansion  of  the  interaction  potential,  with  the  remainder,  results 


(f/K  |v, 


q  (  in  ' 

me  *'K2e0nrcocJ 


1/2 

(^/(z)|e-p|^.(z)). 


(2.58) 


Given  that  the  envelope  function  only  depends  on  z,  only  the  z  component  of  the  dot 
product  ep  applies  to  resulting 


qti  (  Ih  j1' 
me*  (z)y2e0nrcoc  ^ 


cos(^) 


(2.59) 


where  (f)  is  the  angle  between  the  vector  e  and  the  z  axis. 

The  insertion  of  Equation  (2.58)  into  Equation  (2.42)  results  in  the  rate  of 
transition  per  state,  given  by 

2 

w  =  cos \<t>)S{Ef  - Et  -  hco) .  (2.60) 

(m*y  eQnrcoc  \  dz  / 

Notice  that  the  delta  function  in  this  case  does  not  depend  on  K  (the  subbands  are  consid¬ 
ered  identical  parabolas  displaced  by  Ef  -  Ej ). 

At  this  point  it  is  necessary  to  consider  all  occupied  states  in  the  initial 
subband  and  the  correspondent  available  states  in  the  final  subband,  introducing  the  oc¬ 
cupation  statistics  and  integrating  over  all  states.  The  result  is  given  by  the  total  number 
of  carriers,  the  appropriate  occupation  statistics,  and  the  oscillator  strength.  Mathemati¬ 
cally, 

w  =  ,’tqL  It>2  (<r A*) $-<?,(*))  CO s\mEf-E-t>a>) 

(me*)  £0nrcoc  \  dz  /  (2.61) 

x  ALdf™  (Ei)[l-ffFD(Ef)']. 
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Recall  that  d  represents  the  doping  density.  For  low  temperature  operation,  the  total  num¬ 
ber  of  carriers  can  be  considered  d  times  the  volume  AL.  The  oscillator  strength  is  com¬ 
monly  used  as  a  measure  of  the  strength  of  the  transition  and  is  defined  as  [27] 

2 

cos2(0),  (2.62) 
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y/Az) 
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and  obeys  the  rule  [27] 

IX  =1-  (2-63) 

/ 

This  condition  is  valid  for  all  initial  states  i  and  the  sum  of  the  extending  over  all  final 
states  f. 

Finally,  the  absorption  coefficient,  for  transitions  between  two  bound 
states  within  the  conduction  band,  can  be  evaluated  through 

£V.(S‘">  =  -  E,  -  trn) .  (2.64) 


Ideally  the  bound-to-bound  intersubband  absorption  occurs  in  a  single 
wavelength;  however,  in  practice,  the  absorption  spectrum  is  modified  by  various  factors 
[26]  such  as  nonuniformity  in  well  width  and  collisions  experienced  by  electrons,  among 
others.  The  effect  of  this  broadening  is  represented  by  replacing  the  delta  function  by  a 
Lorentzian  function  given  by 


g(hco)  = 


J _ 

In  (Ef 


_ T _ 

Ei-hco)2  +  { r/2)2  ’ 


(2.65) 


where  T  is  the  broadening  parameter,  which  is  equal  the  full  width  at  half  maximum  of 
the  absorption.  The  broadening  parameter  is  normally  determined  experimentally  and  is 
about  10  to  20  meV. 


Further  simplification  is  possible,  when  the  system  is  to  operate  in  low 
temperatures  or  the  difference  between  the  initial  and  final  states  is  much  greater  than 
kBT .  In  this  case,  the  occupancies  probabilities  can  be  considered  one  for  the  initial  state 
and  zero  for  the  final  state,  resulting  in 
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(2.66) 


(Ef-Ei-hco)2  +  (T/2) 


Analyzing  the  oscillator  strength  it  is  possible  to  observe  the  existence  of 
some  selection  rules.  The  matrix  element  shows  that  for  symmetric  wells,  transitions  be¬ 
tween  states  with  the  same  parity  are  not  allowed.  This  can  be  overcome  by  breaking  the 
symmetry  similarly  to  the  interband  case.  Also,  for  normal  incidence  (TE  mode),  cos(^) 
is  zero,  and  consequently,  there  is  no  absorption.  Some  optical  techniques  can  be  applied 
to  overcome  this  selection  rule,  such  as  coupling  throughout  beveled  edges  [28]  using  dif¬ 
ferent  angles  of  incidence,  application  of  several  configuration  of  diffraction  gratings  [29] 
[30]  to  allow  normal  incidence,  introduction  of  strain  in  the  heterostructure  [31],  etc. 

A  strong  characteristic  of  this  kind  of  absorption  is  its  narrow  spectrum 
(about  1  pm).  When  a  broader  absorption  spectrum  is  desired,  it  is  worth  to  consider 
transitions  between  bound  states  to  the  continuum  [32], 


c.  Bound-to-Continuum  Transitions 

Bound-to-continuum  transitions  are  relatively  weak,  when  bound-to- 
bound  transitions  with  high  oscillator  strength  happen  in  the  structure.  This  can  be  ex¬ 
plained  by  Equation  (2.63).  It  becomes  considerable,  however,  when  only  one  confined 
state  exists  in  the  conduction  band. 

To  be  able  to  model  those  transitions,  the  continuum  states  must  be  com¬ 
puted.  A  simple  but  powerful  technique  to  do  that  is  to  introduce  a  fictitious  square  well 
of  width  LF,  and  infinite  barrier  height  within  which  the  continuum  electrons  are  trapped. 
This  configuration  is  shown  schematically  in  Figure  4  [5]. 
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Figure  4  Schematic  diagram  of  the  pseudo-quantification  technique  of  the  continuous 

states  in  of  a  single  square  well. 

The  eigenenergies  of  such  structure  is  given  by  [5] 

,  fi2k2L 

En-n2—^,  (2.67) 

^ me,b 

where  me  b  *  is  the  effective  mass  of  the  electron  in  the  barrier,  and  kLp  is  the  wavevector 
represented  by 

<2-68) 

The  final  states  wavefunctions  are  then,  numerically  calculated  using  the  pseudo¬ 
quantization  model  for  each  En. 

For  example,  if  Lp  is  of  the  order  of  1  pm ,  for  example,  the  energy  level 
separation  close  to  the  barrier  height  is  of  the  order  of  a  hundredth  of  meV,  which  can  be 
considered  much  smaller  than  the  typical  interactions  or  thermal  energies  (order  of  meV). 
The  energies  are  so  close  that,  rather  than  attempting  to  take  into  account  individually, 
they  are  grouped  together  by  means  of  infinitesimal  batches  of  the  one-dimensional  den¬ 
sity  of  states,  given  by  [5] 
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(2.69) 


pw(E)  = 


Lf 

2  nh\  E 


Note  that  as  Lf  tends  to  infinity,  the  density  of  states  increases  without  bound.  More 
states  become  available  over  the  same  range  of  energy  because  the  separation  between 
levels  decreases.  Also,  to  be  consistent  with  the  adopted  energy  reference,  the  barrier 
height  is  subtracted  from  the  final  state  energy,  leading  to 


pw(E)  = 
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(2.70) 


The  density  of  the  continuum  states  must  be  included  in  the  transition  rate 
formulation.  This  procedure  allows  writing  the  expressions  for  the  absorption  coeffi¬ 
cients,  for  both  discussed  cases,  as: 
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(2.72) 


Notice  that  in  this  situation  the  incident  light  energy  can  have  values  that  obey 

tuo  =  Ef-Ei,  for  Ef>V0.  (2.73) 

Also,  the  dependence  of  the  absorption  coefficients  on  LF  is  removed,  since  the  contin¬ 
uum  states  wavefunctions  are  normalized  by  a  factor  proportional  to  the  square  root  of 
the  quantity  LF.  This  fact  is  very  important  because  it  assures  that  the  computation  of  the 
bound-to-continuum  absorption  coefficients  does  not  depend  on  the  model. 
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7.  Addition  of  an  Electric  Field 

Actual  quantum  well  photodetectors  operate  in  the  photoconductive  mode,  de¬ 
manding  the  application  of  bias  through  the  structure  to  be  able  to  extract  the  photocur¬ 
rent.  Under  this  condition,  a  new  potential  factor  is  added  to  the  Hamiltonian  described  in 
Equation  (2.1),  that  becomes  [26] 

H  =  +  Vc(r)  +  Ve(r)  -  qFz ,  (2.74) 

2  m0 

where  the  last  tern  represents  the  potential  change  caused  by  the  applied  electric  field  F. 
Notice  that  if  F  is  constant,  the  overall  potential  bends  linearly,  as  depicted  in  the  Figure 
5. 


Figure  5  Potential  profile  diagrams  of  symmetric  and  asymmetric  quantum  wells  under 

the  effect  of  a  constant  electric  field  F  (bias). 

The  solutions  to  the  Schrodinger’s  equation,  considering  the  new  Hamiltonian  gives  the 
envelope  wavefunction  [18] 

yf(z')  =  C,  Ai(z')  +  C2  Bi(z') ,  (2.75) 

where  Ai  and  Bi  are  Airy  functions,  C;  and  C2  are  constants  that  can  be  evaluated  using 
the  boundary  conditions,  and  z  is  given  by  [18] 
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It  is  important  to  mention  that  the  Airy  functions  present  some  difficulties  in  being  evalu¬ 
ated  numerically  for  either  very  large  or  very  small  arguments. 

The  effect  of  the  static  electric  field  in  the  confined  energy  levels  within  a  quan¬ 
tum  well  is  known  as  quantum  confined  Stark  effect  and  can  be  evaluated  using  the  time 
independent  perturbation  theory  [17]. 

To  determine  the  change  in  the  difference  between  the  fist  and  the  second  con¬ 
fined  energy  levels,  under  the  influence  of  low  electric  field  (case  of  interest  in  this 
work),  a  dipole  perturbation  Hamiltonian  of  the  form, 

Vb(z)  =  -qFz,  (2.77) 


is  employed. 


When  the  quantum  well  is  symmetric,  the  perturbation  is  null  to  the  first  order. 
The  second-order  effect  is  given  by  [5] 


A E  =  2  q2F2 


E\-e2  ’ 


(2.78) 


where  A E  is  always  positive  indicating  a  blue  shift  in  the  photon  absorption.  This  effect 
is  very  small  and  is  normally  obscured  by  the  linewidth  of  the  absorption  spectrum. 

In  asymmetric  quantum  wells  the  first-order  perturbation  is  non-zero  and  ex¬ 
pressed  by 

A E  =  qF[(^  \z\y/x)-(ii/2  \z\v2)\ ,  (2.79) 


where  the  term  between  rectangular  brackets  represents  the  displacement  of  the  average 
position  of  the  electron  from  the  first  to  the  second  confined  state.  In  this  case,  A E  can 
be  either  positive  or  negative,  depending  on  the  bias  signal  and  the  well  geometry.  For 
the  well  depicted  in  Figure  5,  a  blue  shift  occurs  when  positive  bias  is  applied  and  a  red 
shift  occurs  for  negative  bias.  In  asymmetric  structures  the  Stark  effect  becomes  signifi¬ 
cant  and  it  is  commonly  used  in  tunable  detectors  [33], 
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Perturbation  theory  breaks  down  when  the  electron  wavefunction  extends  away 
from  the  well  due  to  field  assisted  tunneling  and  thus  the  prediction  of  absorption  peak 
position  becomes  more  and  more  inaccurate.  [5] 

The  understanding  of  the  quantum  effects  present  in  heterostructures  formed  by 
periodic  repetition  of  ultra-thin  layers  of  semiconductors  with  different  bandgaps  is  nec¬ 
essary  step  in  order  to  design  QWIP  to  meet  specific  requirements.  Also  it  is  important  to 
understand  how  these  heterostructures  can  act  as  photodetectors.  This  section  gathered  a 
set  of  theoretical  tools  that  facilitates  the  comprehension  of  the  discussed  semiconductor 
configuration  in  a  general  sense.  The  following  discussion,  on  the  other  hand,  explains 
how  QWIPs  work  and  compiles  the  mathematical  models  to  evaluate  and  analyze  their 
performance. 

B.  QUANTUM  WELL  INFRARED  PHOTODETECTORS  (QWIP) 

Quantum  well  photodetectors  are  basically  periodic  repetitions  of  the  bar¬ 
rier/well/barrier  structure  between  heavy  doped  contact  layers  grown  epitaxialy  on  a 
semiinsulating  substrate.  Ideally  it  is  possible  to  pile  as  many  stacks  as  the  desired  detec¬ 
tion  bands.  In  practice  up  to  four  bands  have  been  demonstrated  [34],  In  this  work,  a  de¬ 
tector  consisting  of  three  stacks  with  different,  uncoupled,  periodic  repeated  well/barrier 
structure,  separated  by  heavy  doped  contact  layers  is  considered. 

When  a  quantum  well  photodetector  is  exposed  to  an  incident  photon  flux,  transi¬ 
tions  between  quantized  energy  levels  can  occur,  as  mentioned  previously.  The  photon 
flux  can  be  sensed  when  bias  is  applied  to  the  structure  and  the  transitions  allow  the  ex¬ 
traction  of  the  electrons  from  the  well,  to  the  barrier  regions  leading  to  a  photocurrent. 
When  the  final  state  is  confined,  the  electrons  must  tunnel  through  the  potential  barriers, 
made  thinner  by  the  applied  electric  field.  Once  the  electron  has  escaped  the  confines  of 
well,  it  accelerates  toward  the  positive  contact  until  it  is  captured.  Figure  6  shows  a 
schematic  diagram  of  a  biased  multi  quantum  well  structure  and  the  photocurrent  genera¬ 
tion  process. 
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Figure  6  Schematic  diagram  of  the  photocurrent  generation  process  in  a  biased  multi 

quantum  well  structure. 

If  the  exited  state  is  near  the  potential  barrier,  which  is  true  in  most  practical 
cases,  the  probability  of  an  electron  to  tunnel  through  the  barrier,  narrowed  by  the  effect 
of  the  bias  voltage,  can  be  considered  nearly  unity  and  independent  of  the  well  position 
[35].  Thus,  a  single  model  is  used  to  estimate  the  photocurrent  for  both  cases,  bound-to- 
bound  and  bound-to-continuum  transitions. 

It  is  assumed  that  the  absorption  occurs  only  in  the  wells  and  the  light  gets  in  the 
structure  (TV  wells),  crosses  the  nth  well,  reaches  the  top  metal  layer  and  reflects  back, 
crossing  the  same  well  twice.  Then,  if  the  structure  is  under  an  applied  electric  field,  F, 
the  photocurrent  can  be  written  as  [36] 

N  wZ 

1 p(F)  =  +  e~a(2N~n)Lw)ale  V(F)T .  (2.80) 

Here  a  is  the  absorption  coefficient,  <f>0  is  the  incident  optical  power,  ha)  is  the  photon 
energy,  L  is  the  period  of  the  multiple  quantum  wells,  Lw  is  the  width  of  the  well,  v(F) 
is  the  drift  velocity  under  an  applied  electric  field  F,  and  r  is  the  exited  carrier  lifetime. 
Since,  typically  [36],  for  quantum  well  detectors,  aLw  «:  1 ,  Equation  (2.80)  can  be  sim¬ 
plified  to 
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Finally,  in  the  case  of  electrons,  the  drift  velocity  is  given  by  [19] 

v(F )  =  ^"F 
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(2.82) 


where  /un  is  the  electron  mobility  and  vs  is  the  electron  saturation  velocity.  Both  parame¬ 
ters  can  be  found  in  [19]  for  the  main  binaries  and  its  ternaries. 

At  this  point  it  is  worth  to  mention  that  under  high  electric  field,  the  kinetic  en¬ 
ergy  of  the  first  electron,  in  the  current  generation  process,  can  be  large  enough  to  cause 
the  transition  of  a  second  electron,  after  collision.  This  process,  repeated  many  times 
down  the  chain  of  wells,  causes  avalanche  multiplication,  increasing  the  photocurrent. 
This  phenomenon  is  detailed  in  [35]  and  [37]. 

The  way  to  characterize  and  measure  the  performance  of  a  photodetector  is 
throughout  its  figures  of  merit.  The  most  common  used  figures  of  merit  for  quantum  well 
detectors  are  discussed  in  the  following  section. 


C.  FIGURES  OF  MERIT 

1.  Absorption  Spectra 

The  absorption  coefficient  is  the  key  parameter  used  in  photodetector  design.  The 
crystal  heterostructures  are  optimized  for  the  absorption  spectra  requirements  and  the  ab¬ 
sorption  coefficient  is  normally  measured  after  the  sample  growth,  before  the  device  fab¬ 
rication. 

A  detailed  discussion  is  presented  in  Section  II.A.6,  where  the  absorption  coeffi¬ 
cients  for  the  cases  of  interest  in  this  work  are  derived.  It  is  worth  to  notice  though,  that 
the  absorption  coefficient  for  quantum  wells  with  bound-to-continuum  transitions  is  sig¬ 
nificantly  lower  and  wider  than  for  those  with  bound-to-bound  transitions.  This  fact  can 
be  attributed  to  the  conservation  of  the  integrated  oscillator  strength.  Mathematically, 
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(2.83) 


a„ 


A  A/A 


constant , 


where  ap  is  the  absorption  coefficient  peak,  A/1//1  is  the  absorption  bandwidth  and  d  is 
the  well  doping  density. 

Experimental  measurement  of  the  absorptance  using  Fourier  Transform  Infrared 
spectrometry  permits  to  obtain  the  absorption  spectra.  This  procedure  is  deeply  explored 
in  Chapter  IV. 


2.  Quantum  Efficiency  and  Photoconductive  Gain 

The  photocurrent  in  photoconductors  is  defined  as  [37]: 


I  =qV-±G, 
hco 


(2.84) 


where  77  is  the  quantum  efficiency  and  G  is  the  photoconductive  gain.  Comparing  Equa¬ 
tions  (2.84)  and  (2.81),  it  is  easy  to  identify  that 

77  =  2  al  (2.85) 


for  al  cl,  and 

N  iL 

G  =  Yjev(F)T.  (2.86) 

72=1 


For  moderate  electric  fields,  which  is  the  case  in  most  practical  cases,  the  mean 
free  path,  v(F)t  ,  is  considerably  larger  than  the  period  of  the  multiple  quantum  wells  (L) 
and  the  summation  of  Equation  (2.86)  can  be  evaluated  as 


L/v(F)  ’ 


(2.87) 


that  represents  the  rate  between  the  exited  carrier  lifetime  and  the  transit  time. 
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3.  Responsivity 

The  responsivity  figure  quantifies  the  amount  of  photocurrent  generated  per  watt 
of  incident  radiant  photon  power.  It  can  be  expressed  mathematically  as 


R(F)  = 


(2.88) 


Combining  the  expressions  (2.88),  (2.84),  (2.85),  and  (2.87)  results  in 


R(F)  =  (led) 


q  v(F)t 
hco  L 


(2.89) 


It  is  important  to  notice  the  responsivity  dependence  on  the  absorption  coefficient 
and  bias,  through  the  drift  velocity.  The  absorption  coefficient  is  responsible  for  the 
shape  (spectrum)  and  the  electric  field  for  the  amplitude.  Increasing  the  availability  of 
carriers  (doping  concentration),  the  photocurrent  will  increase  and  consequently  the  re¬ 
sponsivity. 


4.  Dark  Current 

A  biased  photodetector,  when  no  light  is  incident,  exhibits  dark  current.  Three 
dark  current  generation  mechanisms  in  quantum  well  devices  can  be  easily  identified 
[38].  First,  sequential  resonant  tunneling  can  happen,  causing  electrons  to  “jump”  from 
well  to  well,  through  the  barriers.  This  process  is  independent  of  temperature  and  is  the 
dominant  source  of  dark  current  at  very  lower  temperatures.  It  can  be  reduced  if  larger 
barriers  are  used.  The  second  mechanism  is  thermally  assisted  tunneling  which  involves 
thermal  excitation  and  tunneling  through  the  tip  of  the  barrier  into  the  transport  states. 
This  process  is  the  dominant  source  at  medium  temperatures  and  can  be  reduced  by  plac¬ 
ing  the  final  state  as  far  as  possible  to  the  initial  state,  reducing  the  probability  of  ther- 
moionc  transitions.  The  third  mechanism  is  classical  thermoionic  emission  and  it  is  the 
dominant  source  at  higher  temperatures.  To  reduce  this  effect,  deeper  wells  must  be  con¬ 
sidered  in  combination  with  reduced  available  carriers  (less  doping).  Obviously  this  will 
reduce  the  photocurrent  as  well,  and  a  compromise  between  both  must  be  established. 
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In  order  to  estimate  the  dark  current,  it  is  necessary  to  estimate  the  effective  num¬ 
ber  of  electrons  which  are  thermally  exited  out  of  the  well  into  the  continuum  transport 
states  as  a  function  o  bias.  This  is  given  by  [39] 


ill  *00 

n*(F)  =  -£- f  fFD(E)T(E,F ) dE 
nh  L  iEi 


(2.90) 


where  the  term  outside  the  integral  is  the  2-D  density  of  states  divided  by  the  multi  quan¬ 
tum  well  period  L,  and  the  Fermi  distribution  is  given  by  Equation  (2.33)  where  E  is  re¬ 
placed  by  ( E  -  £  )  and  £,  represents  the  bound  ground  state.  The  bias-dependent  tunnel¬ 
ing  current  transmission  factor  for  a  single  barrier  is  represented  by  T ( E ,  F )  and  can  be 
described  by  [40] 
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(2.91) 


Recall  that  Vo  represents  the  barrier  height. 

Equation  (2.90)  accounts  for  both  thermoionic  emission  above  the  barrier  height 
and  thermoionically  assisted  tunneling  below  the  barrier  height.  It  can  be  solved  numeri¬ 
cally  in  order  to  estimate  the  dark  current,  expressed  by 

Id(F )  =  qn*  (F)v(F)A  (2.92) 


where  A  is  the  cross  sectional  area  of  the  device  and  all  other  parameters  are  known. 

The  main  mechanism  to  reduce  the  dark  current,  not  mentioned  before,  is  to  re¬ 
duce  the  temperature  of  operation.  To  verify  the  relationship  between  dark  current  and 
temperature  it  is  possible  to  assume  T (E)  =  0  for  E  <V0  and  T (E)  =  1  for  E  >V0  in 

Equation  (2.90),  which  is  a  good  approximation  for  low  biased  systems.  The  dark  current, 
then,  can  be  written  as 

ID(T)  =  qvA^-kBTe  k°T  ,  (2.93) 
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where  Ec  is  the  cutoff  energy  (V0  -  Et )  and  EF  is  quasi-Fermi  energy  which  can  be 

evaluated  using  Equation  (2.34).  Therefore,  the  dark  current  divided  by  the  temperature  is 
related  to  the  temperature  exponentially,  as 


oc  e 


EczEf 

k„T 


(2.94) 


indicating  that  it  decreases  as  the  inverse  of  the  temperature  increases. 


5.  Detectivity 

The  detectivity,  D*,  of  a  detector  is  a  measure  of  the  smallest  photon  flux  that  can 
be  measured,  and  clearly  depends  on  the  noise  associated  with  the  detector.  It  is  mathe¬ 
matically  defined  as 


x  NEP 


(2.95) 


where  A f  is  the  bandwidth  of  the  integration  fdter,  and  NEP  stands  for  noise  equivalent 
power,  that  is  the  noise  power  that  gives  a  unity  signal-to-noise  ratio  (SNR).  The  last  can 
be  estimated  using  Equation  (2.88)  as 

NEP  =  hmifl  (2.96) 

RP(F) 


where  the  dark  current  noise,  im(F ) ,  is  given  by 


iND(F)  =  (4qID(F)GAffN 


(2.97) 


Finally,  making  A/  =  1  Hz  and  combining  Equations  (2.95),  and  (2.96),  the  peak  detec¬ 
tivity  in  [nWHz/w]  can  be  written  as 

R  ( F} 

Da  *  (F)  -  si  (2.98) 

lND 


A  more  relevant  figure  of  merit  than  the  peak  detectivity  is  the  blackbody  detec¬ 
tivity  with  is  obtained  by  calculating  the  spectral  overlap  of  the  responsivity  with  a 
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blackbody  emission  spectrum  at  temperature  TBb ■  Thus,  the  blackbody  detectivity  be¬ 
comes 

— Vd,  (2.99) 

lND 

where  the  bias  dependence  was  omitted  for  convenience,  and  Rbb  is  given  by 

\^R(A)WBB(A)dA 

Rbb=-^ji - •  (2.100) 

[,  WBB(A)dA 

Here  the  power  radiated  per  unit  wavelength  interval  at  wavelength  A ,  per  unit  area  of  a 
blackbody  at  temperature  TBB  (blackbody  spectral  density)  is  given  by 

3.71  ch2  1 

Wbb{A)  =  —a 7  rjj  •  (2.101) 

\^AkBTBB  )  _ ^ 

It  is  worth  noting  that,  for  most  applications,  the  blackbody  responsivity  RBB  is 
reduced  only  a  relatively  small  amount  from  the  peak  value  Rp.  Also,  QWIP  detectivity 
increases  almost  exponentially  with  decreasing  temperature  [41].  This  fact  can  be  ex¬ 
plained  by  the  strong  dependence  of  the  dark  current  on  thermoionic  emission  and  ther¬ 
moionic  assisted  tunneling.  Also  it  must  be  considered  that  the  improvement  on  detectiv¬ 
ity  implies  reduction  in  dark  current  noise. 

Another  important  aspect  is  the  dependence  on  the  detectivity  on  doping.  Apply¬ 
ing  the  same  simplification  used  to  determine  the  dark  current  dependence  on  tempera¬ 
ture,  it  has  been  proved  [28]  that 

D*  °c  .  p  (2.102) 

s]ep  - 1 

where  p  is  linearly  proportional  to  the  doping  density  in  the  well.  This  indicates  that  the 
detectivity  will  increase  with  the  increasing  doping  density  until  reaches  a  maximum  and 
starts  to  decrease  again.  It  is  worth  mentioning  that  the  maximum  region  is  very  wide,  al¬ 
lowing  the  doping  density  to  be  used  to  optimize  the  responsivity  instead. 
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Although  the  mathematical  models  discussed  above  rely  on  a  great  number  of 
simplifications,  it  has  been  shown  (in  the  references  cited  in  this  section)  that  they  are  ap¬ 
propriate  for  the  practical  situations  of  interest. 

The  figures  of  merit  discussed  in  this  section  should  be  considered  in  both  the  de¬ 
sign  and  the  characterization  phases.  In  the  design  they  act  as  requirements  and  great  care 
is  necessary  since  the  maximization  of  one  figure  may  imply  the  minimization  of  other. 

In  the  characterization  phase  they  are  a  measure  of  the  performance  and  the  requirement 
compliance. 


This  chapter  assembled  the  minimum  theoretical  background  necessary  to  build 
the  mathematical  models  that  allows  us  to  design  and  analyze  the  performance  of  QWIP. 
Specificities  of  the  photodetectors  in  discussion  were  considered  and  used  as  simplifica¬ 
tion  factors  in  the  models.  Supported  by  this  theory  the  next  chapter  presents  the  practical 
considerations  towards  the  actual  devices,  the  computational  tools  developed  to  support 
the  design  and  the  final  proof-of-concept  QWIP  design  that  meets  the  requirements  stated 
in  the  introductory  chapter. 
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III.  QWIP  DESIGN 


The  previous  chapter  covered  the  theoretical  aspects  of  quantum  well  photodetec¬ 
tors.  Nevertheless,  to  be  able  to  design  actual  devices  it  is  necessary  to  perform  a  quanti¬ 
tative  analysis  of  the  physical  behavior,  already  described.  This  is  done  sequentially  in 
this  chapter,  introducing  numerical  solutions  of  the  mathematical  models,  adapted  for  the 
design  specificities  of  the  detectors  requirements. 

The  first  requirement  for  a  photodetector  is  the  operation  spectrum.  As  specified 
before,  the  interest  here  is  to  have  a  QWIP  capable  to  detect  simultaneously  wavelengths 
within  three  distinct  bands:  0.9  to  1.4  pm  (near  infrared);  3.8  to  5.0  pm  (mid  infrared); 
and  8.0  to  10.0  pm  (long  wavelength  infrared). 


A.  INITIAL  CONSIDERATIONS 

There  is  no  established  sequence  in  QWIP  design.  Several  analyses  and  consid¬ 
erations  must  be  done  and  most  are  interactive,  forcing  the  designer  to  go  back  and  forth. 
This  section  discusses  the  preliminary  assumptions  and  analysis  towards  the  basic 
photodetector  configuration. 

Many  ways  have  been  used  to  achieve  multiband  detection  using  quantum  wells 
[28].  However,  currently,  two  ways  predominate,  coupled  quantum  wells  [42]  and  multi¬ 
stack  quantum  wells  [34],  The  former  has  the  advantage  of  using  a  fixed  bias  voltage  to 
detect  all  bands  simultaneously.  Conversely  it  is  hard  to  distinguish  those  without  further 
spectral  analysis.  The  requirement  of  discriminating  between  the  bands  demands  sepa¬ 
rated  readouts,  possible  to  implement  in  the  last  configuration. 

The  multi-stack  of  uncoupled  wells  configuration,  selected  in  this  work,  is  shown 
schematically  in  the  Figure  7.  Three  different  stacks  of  quantum  wells  with  different  ab¬ 
sorption  characteristics  are  piled  between  contact  layers.  The  readout  bias  can  be  opti¬ 
mized  for  each  band  independently.  It  is  important  to  mention,  though,  that  this  configu¬ 
ration  demands  high  accuracy  of  epitaxial  growth  due  to  the  large  number  of  different 
layers  needed. 
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Figure  7  Schematic  diagram  of  a  multi-stack  quantum  well  photodetector  with  inde¬ 
pendent  readouts. 

Considering  this  configuration,  to  minimize  the  dark  current  and  maximize  re- 
sponsivity  and  detectivity,  the  second  bound  state  in  the  conduction  band  (final  state)  is 
more  likely  to  be  placed  just  bellow  the  top  of  the  well  ( Ef  ~  V0 ). 

The  first  well  stack  is  responsible  for  detecting  the  near  infrared.  The  lowest 
wavelength  detected  due  to  intersubband  transitions  reported  on  the  literature  [28]  is  2.7 
pm  .  Thus,  inerband  transitions  must  be  used  for  the  NIR  band.  Furthermore,  according 
to  the  selections  rules  for  interband  transitions,  the  well  must  be  asymmetric  to  allow  the 
transition  the  ground  state  of  the  valence  band  to  the  first  exited  state  of  the  conduction 
band.  Thus,  a  step  well  configuration  seems  to  be  more  apporpriate.  Moreover,  only  in¬ 
trinsic  carriers  participate  in  the  transitions  from  the  valence  to  the  conduction  band, 
eliminating  the  necessity  of  doping  these  wells.  The  main  constraints  here  are  to  control 
the  depth  of  the  well  to  achieve  the  right  wavelength,  as  well  as  to  adjust  the  asymmetry 
to  maximize  the  transition  strength. 

The  second  well  stack  is  responsible  for  mid-infrared  detection.  This  can  be  ac¬ 
complished  by  intersubband  transitions  in  a  deep  quantum  well.  The  main  constraint  here 
is  to  achieve  a  deep  well  without  having  too  much  strain  in  the  system.  Also  the  doping 
concentration  of  these  wells  must  be  adjusted  to  maximize  the  figures  of  merit. 
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The  third  stack  is  responsible  for  the  long-infrared  detection.  This  is  the  most 
commonly  used  configuration  and  can  be  achieved  without  much  dificulty  using  shallow 
wells.  The  principal  constraint  here  is  to  control  the  dark  current.  In  shallow  wells,  the 
thermoionic  emission  is  strong;  hence  the  doping  concentration  must  be  reduced  to  con¬ 
trol  this  process. 

All  the  three  stacks  are  schematically  represented  in  the  Figure  8,  in  terms  of  po¬ 
tential  diagram.  It  is  worth  noting  that  all  considerations  from  now  on  focus  on  the  opti¬ 
mization  of  the  performance  of  the  three  stacks  together.  Since  this  work  aims  a  verifica¬ 
tion  of  the  concept,  a  few  compromises  on  the  requirements  are  acceptable  to  make  it  fea¬ 
sible. 
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Figure  8  Schematic  diagram  of  the  potential  profile  of  the  basic  detector  configuration, 

including  all  three  quantum  well  stacks. 

In  selecting  the  configuration,  the  materials  must  be  chosen,  based  upon  their 
properties,  availability,  growing  process  maturity,  and  cost,  among  other  factors.  Figure  9 
shows  a  plot  of  energy  bandgaps  vs.  lattice  constants  for  major  III-V  compounds  [43], 
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Ternary  alloys  are  denoted  by  lines  between  binaries  and  quaternary  alloys  are  repre¬ 
sented  by  the  colored  areas. 


Figure  9  Plot  of  energy  bnadgaps  vs.  lattice  constants  for  major  III-V  compounds.  The 
solid  lines  and  the  dashed  lines  represent  the  direct-gap  and  indirect-gap  re¬ 
gions,  respectively.  Ternary  alloys  are  denoted  by  lines  between  binaries  and 
quaternary  alloys  are  represented  by  colored  areas.  (After  Ref.  [43].) 

The  requirement  of  closely  matching  lattice  constants  separates  these  semicon¬ 
ductors  into  five  distinct  families  in  terms  of  compatible  substrate  materials  [5],  The  best 
controlled  of  these  families  is  GaAs/AlxGai_xAs,  and  all  the  compositional  aluminum 
fractions  are  accessible,  as  the  lattice  constant  does  not  vary  significantly  from  GaAs  to 
AlAs.  GavIn/_xAs  can  be  grown  on  GaAs  substrates  if  the  composition  of  indium  is  not 
very  high  ( x  <  0.4 ),  and  the  strain  effects  for  this  system  are  well  known  (see  for  exam¬ 
ple  [44]).  Another  important  system  includes  IntGay^AsyP/.j,  and  AlxIni_xAs,  which  can 
be  grown  on  InP  substrates  and  is  commonly  used  in  lasers  applications  [45].  Also, 
InvAs/_vSb  and  AlvGa/_vSb  are  normally  grown  on  GaSb  substrates.  Finally,  the  lattice 
constants  available  to  Flg/_ACdvTe  are  suitable  to  deposit  on  a  CdTe  substrate  and  allow  a 
large  range  of  accessible  bandgaps.  Recently  more  combinations  have  been  tested  [28], 
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In  addition  to  the  lattice  constant  and  bandgap  energies,  the  band  offsets  are  ut¬ 
terly  important  to  determine  the  band  edge  potentials  and  consequently  the  depth  of  the 
wells.  Some  theoretical  models  are  available  to  roughly  estimate  the  band  offsets;  al¬ 
though,  to  design  quantum  well  devices,  much  more  accurate  values  than  those  obtained 
from  the  theoretical  models  are  required.  Methods  such  as  absorption  measurement,  pho¬ 
toluminescence  measurement,  and  x-ray  core  level  photoemission  spectroscopy  are  used 
for  experimental  determination  of  band-offsets  [26],  Even  experimental  methods  cannot 
assure  accurate  values  of  band  offset,  making  this  parameter  a  considerable  constraint 
when  working  with  certain  families  of  materials  [45].  A  good  collection  of  band-offset 
parameters  can  be  found  in  [20]. 

The  analysis  on  available  data  from  the  families  of  compounds  mentioned  above, 
including  cost  and  fabrication  potentialities  led  to  two  possibilities,  GaAs/AlxGa;.xAs/ 
lnvGa/-.vAs  on  GaAs  substrate  and  ( I n o  j?G ao, ^ 7 A s ) ,,( In P ) / .J  Alo.45Ino.52As  on  InP  substrate. 
These  families  are  highlighted  in  red  in  Figure  9  For  simplicity,  from  now  on,  both  are 
addressed  by  their  substrate,  as  the  GaAs  system  and  the  InP  system. 

Very  reliable  material  data  (properties)  and  mature  and  cheaper  fabrication  proc¬ 
esses  are  advantages  of  the  GaAs  system.  However,  to  assemble  deep  wells  necessary  for 
mid-infrared  detection  it  is  necessary  to  increase  the  Indium  composition,  raising  the 
strain  in  the  system  due  to  the  lattice  mismatch.  Limitations  on  indium  composition  are 
well  documented  in  [44]  [46]  and  [47],  and  lead  to  a  design  limit  of  0.25  if  used  in  thin 
wells.  Furthermore,  the  critical  layer  thickness  that  represents  the  maximum  thickness  of 
the  GaAIn /_vAs  under  which  the  large  lattice  mismatch  is  totally  accommodated  by  uni¬ 
form  elastic  strain,  is  dependent  on  the  indium  quantity  [45].  Therefore,  the  higher  the  in¬ 
dium  composition,  the  lower  the  critical  thickness,  limiting  the  design  possibilities.  Also, 
since  AlAs  is  an  indirect  semiconductor,  the  aluminum  concentration  must  be  kept  lower 
than  0.4,  to  have  a  direct  AlvGa/.vAs  ternary.  Moreover  the  high  bandgap  of  GaAs  and  the 
limitation  on  the  indium  composition  indicates  that  it  might  be  difficult  to  comply  with 
the  near  infrared  requirement. 

The  InP  system,  on  the  other  hand,  presents  very  controversial  band  offset  data 
[45],  and  uses  more  difficult  and  expensive  fabrication  processes.  The  absence  of  strain, 
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due  to  the  matched  lattice,  is  an  advantage.  Also  the  lower  bandgaps  of  the  compounds 
indicates  the  possibility  to  achieve  the  desired  near-infrared  detection  spectrum.  Adapta¬ 
tions  on  the  requirements  are  also  considered  in  order  to  make  it  feasible. 

A  significant  limitation  concerning  the  material  composition  is  the  growing  proc¬ 
ess.  When  using  MBE,  for  example,  a  limited  amount  of  the  basic  compounds  vapor 
sources  are  available.  Therefore  the  number  of  different  compositions  of  the  same  ternary 
or  quaternary  is  limited.  For  example,  using  two  sources  of  A1  in  different  temperatures, 
it  is  possible  to  grow  Al<Ga/_<As  with  3  different  fixed- values  of  x  determined  by:  the 
first  source  (when  only  its  shutter  is  open),  the  second  source  (when  only  its  shutter  is 
open),  and  both  (when  both  shutters  are  open).  This  restricts  the  structure  in  having  only 
3  different  compositions  of  AfGa/_xAs,  limiting  the  designing  degrees  of  freedom. 

Lastly  it  is  necessary  to  consider  the  angle  of  incidence  selection  rules  to  be  able 
to  design  the  optical  coupling.  According  to  the  previous  chapter,  intersubband  transi¬ 
tions  are  forbidden  in  symmetric  quantum  wells,  thus  diffraction  gratings  must  be  etched 
on  the  top  of  the  last  contact  layer,  just  before  the  metallization.  Also  a  high  purity  non¬ 
alloy  ed  gold  and  silver  must  be  considered,  in  order  to  maximize  the  reflection  [48],  A 
complicating  factor  is  that  the  grating  must  be  optimized  for  two  different  bands.  Several 
configurations  have  been  proposed  [49]  [50],  however  the  limitations  on  mask  lithogra¬ 
phy  and  processing  capabilities  restrict  the  possibilities. 

Another  complicating  factor  is  the  fact  that  the  wavelengths  of  the  LWIR  are  al¬ 
most  two  times  the  wavelengths  of  the  MWIR.  As  will  become  clear  in  the  Section  IV.2, 
it  is  impossible  to  optimize,  simultaneously,  the  grating  pattern  for  wavelengths  that  are 
multiples  to  each  other.  Therefore,  a  good  strategy  is  to  place  the  detection  peaks  closer 
of  each  other.  Then,  considering  that  smaller  dimensions  are  required  for  smaller  wave¬ 
lengths  and  that  the  small  dimensions  (less  than  1  pm )  are  difficult  to  implement,  the  de¬ 
tection  peaks  for  the  MIWR  and  LWIR  should  be  placed  as  close  to  each  other  as  possi¬ 
ble.  This  is  also  an  advantage  if  just  one  grating  pattern  is  possible  to  be  applied  to  satisfy 
both  spectral  regions. 


44 


Having  established  the  basic  detector  configuration,  based  on  the  requirements, 
constraints,  and  data  availability,  the  first  design  step  is  to  quantitatively  determine  the 
bound  states  in  the  conduction  and  valence  bands  and  their  respective  wavefunctions. 


B.  SELF  CONSISTENT  SCHRODINGER-POISSON  SOLUTIONS 

The  self-consistent  Schrodinger-Poisson  solutions  are  necessary  to  account  on  the 
potential  distortion  caused  by  the  charge  distribution.  This  approach  can  be  summarized 
as  follows.  First,  the  band  edge  potential  Ve(z)  is  computed  from  the  structure  material 

properties.  Then  the  perturbation  potential  Vb(z) ,  due  to  applied  electric  field,  is  com¬ 
puted  and  added  to  the  band  edge  potential  (V(z)  =  Ve(z)  +  Vh(z) ).  Next  step  is  to  solve 
the  Schrodinger’s  equation  for  the  structure,  to  obtain  the  bound  states  Et  and  their  re¬ 
spective  wavefunctions  y/i ,  that  will  be  used  as  the  first  guess  to  compute  the  charge  dis¬ 
tribution  potential  Vp{z ) ,  through  the  solution  of  the  Poisson’s  equation  of  the  structure. 
The  potential  V p(z)  is  then  added  to  the  original  (V(z)  =  Ve(z)  +  Vh(z)  +  V p )  and  the 

Schrodinger’s  equation  is  solved  to  the  new  potential.  The  iteractions  repeat  until  the 
bound  states  energies  converge.  This  procedure  is  schematically  depicted  in  Figure  10. 
The  details  in  the  figure,  including  equations  and  the  respective  parameters  are  discussed 
in  the  following  sections. 
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Figure  10  Schematic  diagram  of  the  self-consistent  Schrodinger-Poisson  solution  main 
steps.  All  equations  and  parameters  are  explained  on  the  text. 


The  first  step  to  solve  a  quantum  well  Shrodinger’s  equation  is  to  compute  the 
band  edge  potential  profile.  This  can  be  done,  basically,  taking  the  bandgap  energies  of 
the  material  layers  along  with  their  band  offset  ratios. 


1.  Bandgap  and  Effective  Mass  Parameters 

The  bandgap  parameters,  for  both  the  GaAs  system  and  InP  system,  shown  in  the 
Figure  9  are  tabulated  in  [20].  For  the  ternary  AXB{_XC ,  the  following  expression  is  used: 


Eg(AxB,_xC )  =  xEg(AC)  +  (1  -  x)Eg(BC )  -  x(l  -  x)Cb(ABC )  (3.1) 
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where  the  so  called  bowing  parameter  Cb  accounts  for  the  deviation  from  a  linear  interpo¬ 
lation  between  the  endpoints  two  binaries  AC  and  BC.  For  the  lattice  matched  quaternar¬ 
ies  ( AxBl_xC)y(DE)l_y ,  the  corresponding  band  parameter  is  expressed  as 

Es  ((AJ)l-IQ,(DE)I_y )  =  yE/JxB,_xC)  +  (1  -  y)Eg(DE)  -  y(  1  -  y)Cb(ABCDE) .  (3.2) 

The  temperature  dependence  of  the  bandgap  must  also  be  considered  and  is  given 
by  Equation  (2.13).  Thus,  the  substitution  of  the  GaAs  and  InP  systems  material  parame¬ 
ters  into  Equations  (2.13),  (3.1)  and  (3.2)  gives  the  needed  bandgap  values. 

The  effective  mass  profile  is  directly  deduced  by  the  properties  of  each  layer.  For 
the  binary  compounds,  the  effective  masses  along  the  crystal  direction  [100]  are  given  by 
the  following  expressions  [20]: 


i  =  (i +2Fr)+MVL(^M 

Me*  AJ 


-^  =  *-2  r2,  (3.4) 

mhh* 

and 

=  Yi  +  2f2  (3-5) 

™,h * 

where  FK  is  the  Kane  parameter,  E  is  the  Kane  energy,  Aso  is  the  spin-orbiting  splitting 
parameter,  and  yi  are  the  Luttinger  parameters.  All  these  parameters  are  obtained  ex¬ 
perimentally  and  are  listed,  for  the  binaries  in  question,  in  [20],  For  the  ternaries  and  lat¬ 
tice  matched  quaternaries  of  the  type  A  VB  i_x,  the  effective  masses  can  be  obtained  by  the 
linear  interpolation  of  the  extremes  (A  and  B)  through  the  expression 

- 1 - =  x — - —  +  (1  —  x) — - — .  (3.6) 

m*(AxBl_x)  m*(A)  m*(B) 

The  substitution  of  the  GaAs  and  InP  system  material  parameters  into  Equations  (3.3)  to 
(3.6)  gives  a  series  of  equations  necessary  to  obtain  the  needed  effective  masses. 
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2.  Potential  Profile 

An  important  remark  is  that,  since  the  design  seeks  a  detector  with  uncoupled 
wells,  all  development  in  this  section  considers  a  configuration  of  two  external  barriers 
and  one  quantum  well.  Also,  step  wells  with  one  or  more  steps  can  be  used  to  create 
asymmetric  structures.  In  this  sense,  the  bottom  of  the  conduction  band  of  the  smallest 
bandgap  layer  is  considered  the  potential  zero,  and  the  center  of  the  same  layer  is  taken 
as  z  =  0 .  Thus  the  conduction  band  edge  potential  profile  is  given  by 


KC*(z)  = 


N-l 


£  [(£„.,  -  ■ E„  )<C  ]  @(-  -  z. )  Z  >  o 

n= 1 
-M+ 1 

y[(E  .-E  )SCB  ,l©(-z  +  z  )  for  z  <  0, 

/  j  |_V  g,m- 1  g,m  s  m,m- 1 J  V  m  /  ’ 


(3.7) 


m=—  1 


where  N  is  the  number  of  layers  in  the  positive  direction,  SFB+l  represents  the  conduction 

band  offset  ratio  between  the  nth  and  (n+1  )th  layers,  found  in  [20]  and  [45],  and  zn  is  the 

coordinate  of  the  interface  between  them.  The  same  interpretation  is  taken  for  the  nega¬ 
tive  direction,  replacing  A  by  M  and  n  by  m.  The  bandedge  potential  defined  this  way  can 
handle  several  distinct  layers  to  model  asymmetric  step  wells.  In  the  same  way,  the  va¬ 
lence  band  potential  is  given  by 


v;*(z)  = 


N-l 


-E 


Z  E„  )(i  -  <0]  e(z  -  z„ ) 

n= 1 
-M+ 1 

V  \-{E  ,-E  )(\-SCB  ,)le(-z  +  z  ) 

k  g,m— 1  g,m  /  V  m,m—\  7  J  V  m  ’ 


REF 


for  z  >  0 


m=- 1 


(3.8) 


-EfF  for  z  <  0, 


where  Ebef  is  the  lowest  bandgap  within  the  structure. 

As  mentioned  before,  the  most  critical  parameter  to  determine  the  band  edge  po¬ 
tential  is  the  band  offset.  This  is  determined  experimentally  and,  normally,  gives  values 
widely  spread,  depending  on  the  experimental  method  used,  especially  for  the  InP  sys¬ 
tem.  The  values  used  were  confirmed  experimentally  in  the  Sensor  Research  Laboratory 
(SRL),  Naval  Postgraduate  School  (NPS)  [51],  for  the  GaAs  system,  and  an  average 
value  suggested  by  the  authors  in  [20]  and  [45],  for  the  InP  system.  When  compounds 
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with  slightly  different  compositions  from  the  experimental  data  are  used,  it  is  suggested 
on  the  literature  that  linear  interpolation  would  be  a  reasonable  correction.  The  data  used 
in  this  work  is  given  in  [20], 

Recall  from  Equation  (2.77)  that  the  potential  due  to  bias  applied  on  the  het¬ 
erostructure  is  given  by 

Vb(z)  =  -qFz.  (3.9) 

The  addition  of  this  and  the  band  edge  potential  is  represented  by 

V(z)  =  Vb(z)  +  Ve(z),  (3.10) 

where  V (z)  is  the  potential  profile  without  taking  into  account  the  electrostatic  potential. 

At  this  point,  all  information  is  in  place  to  solve  the  structure  Schrodinger’s  equa¬ 
tion  and  to  obtain  the  confined  energy  levels  as  well  as  the  corresponding  wavefunctions. 


3.  Numerical  Solution  of  the  Schrodinger’s  Equation 

The  selected  method  to  solve  numerically  the  Schrodinger’s  equation  here  was  a 
combination  of  the  shooting  method  and  the  finite  differences  [52]  [18].  The  one- 
dimension  Shrodingers’s  equation  in  the  heterostructure  along  the  growth  axis  is  given  by 

-7T — ^-—^rV(z)  +  V(z)yf(z)  =  Ey/(z),  (3.11) 

2  oz  m* (z)  dz 


that  can  be  expanded  as 


1  d  *  .  d 

- - - m*(z) — y/(z)  + 

[m*(z)]2  dz  dz 


-4-r|W(z)  =  zra-fM')  • 

m  *  (z)  dz  n 


(3.12) 


First,  the  continuous  z  direction  is  discretized  in  an  evenly  spaced  grid  with  step 
Sz .  Second,  the  derivatives  are  replaced  by  their  second  order  central  finite  difference 
approximations,  expressed  as 


df  _  f(z  +  Sz)-f(z-Sz) 
dz  z  2  Sz 


(3.13) 


for  the  first  derivative,  where / represents  an  arbitrary  function,  and 
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(3.14) 


d2/  ..  f(z  +  Sz)  -  2/(z)  +  /(z  -  Sz) 
dz 2  z  (£z)2 

for  the  second  derivative.  Applying  Equations  (3.13)  and  (3.14)  to  Equation  (3.12),  after 
some  algebra,  the  result  is  [18] 


y/{z  +  dz) 
m*  (z  +  Sz  /  2) 


'  2  (Szf 
h2 


[V(z)-E]  + 


1 

m*  (z  +  Sz  /  2) 


y/{z-Sz) 
m*  (z  -  Sz  /  2) 


(3.15) 


Equation  (3.15)  is  the  shooting  equation.  The  equation  shows  that  knowing  two  values  of 
the  wavefunction,  y/(z  -  Sz)  and  y/(z) ,  the  third  point,  y/{z  +  Sz) ,  can  be  predicted.  Us¬ 
ing  the  new  point  together  with  its  predecessor,  a  fourth  point  can  be  calculated  and  so 
on.  Hence  the  wavefunction  can  be  deduced  for  any  particular  energy  value.  The  solu¬ 
tions  for  stationary  states  have  wavefunctions  which  satisfy  the  standard  boundary  condi¬ 
tions,  expressed  as 

0 

^r(z)  — >  0  and  — i//(z)—>  0,  as  z— >±°°.  (3.16) 

oz 

It  has  been  proved  [18]  that  the  general  initial  conditions  that  can  handle  any  po¬ 
tential  profile  is  given  by 

y/{z  -  Sz)  =  0  and  y/(z)  =  1,  (3.17) 

since  multiplying  an  eigenstate  (wavefunction  if/)  for  a  constant  does  not  affect  the  ei¬ 
genvalue  (bounded  energy  E). 

Sequentially,  the  energy  is  varied  systematically  in  SE  steps,  until  the  wavefunc¬ 
tion  switches  from  diverging  to  to  diverging  to  .  That  can  be  mathematically  ex¬ 
pressed  as 

(v{E)\m  ){v(E  +  SE)  1^ )  <  0 .  (3.18) 

Clearly  an  energy  level  exists  between  these  values  for  which  the  wavefunction  will  tend 
smoothly  to  zero.  Since  the  problem  now  is  to  find  a  root  of  \f/{E) |  =  0 ,  the  energy 
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grid  is  refined  in  the  interval  ( E ,  E  +  SE )  and  the  Newton’s  method  is  applied  [52],  This 
is  done  to  accelerate  the  convergence  to  the  true  solution.  Recall  that,  according  to  the 
Newton’s  method,  if  Et  is  the  first  guess  to  the  solution  of  W(.E)\  =  0 ,  then  a  better 

estimate  is  given  by 


V'(E,) L„  ’ 


(3.19) 


where  the  denominator  of  the  quotient  at  the  right  side  is  the  first  derivative  of  the  wave- 
function  f//  with  relation  to  E,  at  Et ,  at  z  — >  . 


The  tolerance  of  1  meV  in  energy  values  is  certainly  within  the  error  caused  by 
other  uncertainties.  In  terms  of  wavelength,  for  a  quantum  well  structure  with  absorption 
peak  around  4.5  pm ,  the  change  of  1  meV  in  the  difference  between  the  conduction  band 
bound  states  would  cause  a  shift  of  the  order  of  one  hundredth  of  a  micrometer  on  the 
peak  absorption. 

In  practice,  the  above  mentioned  iterative  procedure  starts  and  ends  at  finite  dis¬ 
tances  from  z  =  0 .  This  is  the  same  as  assuming  that,  at  the  lower  and  upper  limits  (effec¬ 
tive  infinites)  of  the  z  domain,  the  barrier  potential  is  infinite,  forcing  the  wavefunction  to 
be  zero  there.  The  limits  of  the  z  domain  should  be  chosen  to  be  at  sufficient  extent  so  as 
not  to  affect  the  eigenvalue.  It  has  been  proven  [18]  that  if  the  energy  is  the  only  motiva¬ 
tion  for  this  calculation,  outer  barriers  as  wide  as  150  A  will  suffice. 


When  the  wavefunction  is  needed  to  compute  other  parameters,  great  care  must 
be  taken.  Even  though  the  eigenvalues  can  be  obtained  within  the  convergence  tolerance, 
for  thin  outer  barriers,  the  wavefunctions  might  not  satisfy  the  second  boundary  condi¬ 
tions,  y/'  =  0  at  z  — »  ±oo .  On  the  other  hand,  for  thick  outer  barriers,  the  wavefunction 
starts  to  diverge,  due  to  limitation  in  numerical  precision  of  the  energy  values.  The  defini¬ 
tion  of  thin  and  thick  depend  on  the  structure  being  solved  and  should  be  monitored  for 
each  specific  problem,  in  order  to  find  the  appropriate  limits  of  the  z  domain.  Since  the 
wavefunction  convergence  is  dependent  on  the  accuracy  of  energy,  the  greater  the  figure 
accuracy  that  can  be  used,  the  thicker  the  barriers  that  can  be  included  in  the  model.  Per- 
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sonal  computer  codes  normally  assign  16-figure  accuracy  to  a  double  precision  number. 

If  the  quantum  well  configurations  demand  higher  precision  than  that,  more  sophisticated 
computational  tools  must  be  used. 

The  last  parameter  that  must  be  defined  is  the  step  length  Sz .  Since  the  atomic 
dimension  limits  the  layer  thickness,  there  is  no  practical  reason  to  use  steps  smaller  than 

1A.  Higher  values  can  be  used,  to  reduce  the  computation  time,  according  to  the  configu¬ 
ration  being  solved.  For  the  structures  computed  in  this  work,  Sz=  1A  was  used. 

The  final  wavefunction  is  then  obtained  applying  the  energy  value  calculated  by 
the  Newton’s  method  in  Equation  (3.15),  considering  the  appropriate  limits,  steps  and 
tolerances,  just  discussed. 

The  wavefunctions  obtained  for  this  numerical  method  are  not  normalized.  This 
can  easily  be  achieved  multiplying  the  wavefunction  by  its  normalization  factor, 

«,  =  ,  1  ,  (3.20) 

Jr >*(zmi)4 

V  Jzmin 

where,  the  symbol  (*)  indicates  the  complex  conjugate. 

The  shooting  method,  sketched  in  the  Figure  1 1,  is  a  powerful  tool  to  handle  any 
kind  of  potential  profile.  That  is  very  important  in  solving  biased  multiple  layer  structures 
with  a  varying  doping  profile.  Because  of  that,  it  is  very  suitable  for  self-consistent  solu¬ 
tions.  However,  if  not  carefully  implemented,  it  can  be  very  costly  computationally.  Also 
the  limitation  on  the  wavefunction  convergence  for  wide  outer  barriers  must  be  taken 
care  externally.  It  is  normally  possible,  imposing  a  zero  value  to  the  wavefunction  before 
the  divergence  region. 
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Figure  1 1  Schematic  diagram  of  the  shooting  method  procedure. 


The  shooting  method  was  implemented  in  Matlab  to  obtain  the  confined  energy 
levels  and  its  respective  wavefunctions  of  single-period  heterostructures.  Routines  con¬ 
taining  the  material  properties  of  the  GaAs  and  InP  system  were  also  programmed  to 
support  the  calculations.  To  validate  the  code,  several  different  structures  from  the  litera¬ 
ture  were  solved  and  the  results  were  found  in  good  agreement.  Figure  12  shows  the  con¬ 
duction  and  the  valence  band  confined  energy  levels  and  their  respective  wavefunctions 
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of  a  single  step  well  studied  in  [51].  The  structure  dimensions  and  composition  are  also 
shown  in  the  same  figure.  The  wavefunctions  have  been  rescaled  for  visual  purposes. 


Figure  12  The  calculated  wavefunctions  of  the  confined  energy  levels  of  the  asymmetric 

step  well  studied  in  [51]. 


Figure  13  shows  the  conduction  band  confined  energy  levels  and  their  respective 
wavefunctions  of  a  more  complex  configuration,  studied  in  [53].  The  structure  consists  of 
three  different- size  coupled  quantum  wells  separated  by  thin  barriers.  The  dimensions 
and  composition  are  shown  in  the  same  figure. 
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Figure  13  The  calculated  wavefunctions  of  the  confined  energy  levels  of  the  triple- 

coupled  quantum  well  studied  in  [53]. 
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The  results  show  that  this  implementation  of  the  shooting  method  is  valid  for  the 
configuration  of  interest  of  this  thesis  and,  moreover,  it  can  handle  complex  coupled 
structures.  Also  the  effects  of  the  potential  bending  due  applied  bias  and  charge  distribu¬ 
tion  can  be  included  since  the  method  is  applicable  for  any  potential  profile.  This  is  very 
appropriate  for  self-consistent  solutions  as  shown  in  the  next  section. 


4.  Charge  Distribution  Potential  Evaluation 

The  carrier  distribution  in  a  quantum  well  depends  on  the  quasi-Fermi  level,  EP, 
within  the  well.  The  basic  idea  to  determine  it  was  mathematically  described  in  the  Sec¬ 
tion  II. A.  5.  Since  it  is  desired  that  there  be  only  two  confined  states  in  the  conduction 
band,  using  Equation  (2.34),  it  is  possible  to  express  the  total  number  of  electrons  as  [5] 


A  = 


kj 

nfi1 


f 

EF~E 1 

E p  ~E2 
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*  In 

1  +  e  k>T 

+  /nf  *  In 

1  +  e  k°T 

V 

) 

(3.21) 


were  Ei  and  A  represents  the  first  and  second  bound  states,  respectively,  and  m‘‘v  *  is  the 

average  effective  mass,  accounting  on  the  probability  of  the  electron  location  along  the 
structure,  for  each  energy  level.  It  can  be  evaluated  through  [5] 


!=  1 


*\Vi  *  (z)wXz)  dz 


(3.22) 


where  the  effective  mass  of  the  zth  layer,  nr  * ,  is  weighed  by  the  total  probability  of  find¬ 
ing  an  electron  at  that  layer,  represented  by  the  integral  term.  Since  the  total  number  of 
carriers  can  be  obtained  directly  by  the  doping  concentration,  the  only  unknown  in  Equa¬ 
tion  (3.21)  is  the  quasi-Fermi  energy.  Using  the  iterative  Newton’s  root  finder  to  solve 
the  equation 


<n 

(3q  ^ 

1 

ii 

< 

f 

EF  E] 

Ep—E2 

\ 

m“v*  In 

1  +  e  hT 

+  m“v*  In 

1  +  e  k°T 

V 

- 

- 

) 

where  the  first  guess  of  EF  can  be  obtained  considering,  initially,  only  the  first  state  popu¬ 
lated  ( T  -  0  ),  which  is  given  by  [5] 
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Hence,  a  better  estimate  can  be  obtained  by  repeating  [52] 
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^  F  ^  F 


/(4) 
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(3.24) 


(3.25) 


until  the  convergence.  A  tolerance  of  1  meV  is  suitable  for  this  evaluation. 

Having  calculated  the  quasi-Fermi  level,  the  total  number  of  carrier  populating 
each  subband  is  obtained  using  Equation  (2.34).  For  very  low  temperature  and  deep 
wells,  this  can  be  avoided  and  it  is  reasonable  to  consider  all  electrons  populating  the 
ground  state  (first  subband). 

The  next  steps  are  to  calculate  the  net  surface  charge  density,  the  electric  field 
along  the  structure,  and  finally  the  charge  distribution  potential.  Those  can  be  obtained 
directly  from  Equations  (2.39),  (2.40),  and  (2.36),  respectively. 

The  evaluation  of  the  charge  distribution  potential,  through  the  Poisson’s  equation 
solution,  sketched  in  Figure  14,  becomes  important  in  heavily  doped  structures,  such  as 
the  proposed  in  this  work.  Also,  the  possibility  to  monitor  non-measurable  parameters 
such  as  electric  field  and  charge  density  inside  the  structure,  allows  deeper  physical 
analysis. 
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Figure  14  Schematic  diagram  of  the  charge  distribution  potential  evaluation 


The  charge  distribution  potential  computation  was  included  in  the  Matlab  code  as 
part  of  the  self-consistent  approach.  The  self-consistent  Schrodinger-Poisson  program 
was  validated  through  the  solutions  of  several  structures  from  the  literature.  The  results 
were  found  in  good  agreement.  For  example,  Figure  15  shows  the  potential  profile,  the 
confined  energy  levels  and  their  respective  wavefunctions  for  the  structure  studied  in 
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[54],  Huang  et  al  studied  a  GaAs-based  step  well  with  the  barriers  doped  in  the  middle 
(modulation  doping). 


(a) 


(b) 


Figure  1 5  Self-consistent  potential  profile  with  normalized  wave  functions  for  the  two 
lowest  energy  levels  under  a  30-kV/cm  electric  field  (a)  and  a  -  30-kV/cm 
electric  field  (b),  for  the  structure  studied  in  [54],  Dashed  line  is  the  potential 
profile  without  the  charge  distribution  effects. 


Figure  16  shows  the  individual  potential  contributions  as  well  as  the  doping  profile, 
charge  density  and  electric  field  behaviors  in  the  structure  under  a  30-kV/cm  electric 
field,  corresponding  to  Figure  15  (a). 


Figure  16  Individual  potential  contributions  in  the  structure  (left)  and  the  doping  profile, 

charge  density  and  electric  field  behaviors  (right). 
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At  this  point,  using  the  developed  and  validated  computational  tool,  it  is  possible 
to  calculate  the  confined  energy  levels  and  their  respective  wavefunctions  of  several 
types  of  heterostructures  in  a  wide  variety  of  configurations.  Furthermore,  using  the  ab¬ 
sorption  coefficient  expressions  derived  in  Section  II.A.6,  it  is  also  possible  to  estimate 
the  absorption  coefficient  of  the  quantum  well  systems.  This  was  done  for  both  the  GaAs 
and  the  InP  systems,  considering  the  guidelines  discussed  in  Section  III. A  in  order  to 
maximize  the  absorption  coefficient  on  the  spectral  regions  of  interest.  The  next  section 
brings  this  discussion  and  the  final  designed  structure  configurations. 

E.  DESIGN  OF  THE  QWIP  SEMICONDUCTOR  HETEROSTRUCTURES 

The  computational  tools  developed  to  support  this  work  are  not  optimized  for  de¬ 
sign  purposes.  Optimization  algorithms  must  be  considered  and  more  computational  ca¬ 
pability  must  be  available  to  obtain  that.  However  it  is  possible  to  select  the  configura¬ 
tions  and  material  composition  based  on  the  arguments  listed  in  Section  III.  A  and  to  run 
the  programs  to  analyze  the  results.  Changes  in  the  inputs  must  be  made  based  on  the 
physical  properties  and  the  desired  result.  After  some  interactions  it  is  possible  to  refine 
the  inputs  to  obtain  the  desired  outputs  (requirements).  This  was  done  for  both  systems 
aiming  to  maximize  the  absorption  coefficient,  placing  the  peak  at  the  desired  wave¬ 
length. 


1.  GaAs  System 

In  terms  of  available  technology  and  cost  predictions,  the  GaAs-based  configura¬ 
tion  is  the  most  attractive.  In  view  of  the  guidelines  and  taking  into  account  the  limita¬ 
tions  discussed  in  Section  III.  A,  the  following  design  requirements  were  set  for  the  detec¬ 
tors  based  on  GaAs  system: 

a.  The  peak  absorption  of  the  three  different  bands  for  temperatures  below 

77K  should  be:  the  maximum  value  obtainable  for  NIR  (limited  by  the 
GaAs  bandgap);  as  close  as  possible  to  5.0  pm  for  MWIR  (limited  by 
the  grating  pattern);  and  below  9.0  pm  (limited  by  the  grating  pattern 
and  leakage  current  due  to  thermoionic  emission). 
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b.  Bound-to-quasibound  transitions  should  be  preferred  to  minimize  the 
possibility  of  thermoionic  transitions,  maximize  the  oscillator  strength 
and  reduce  the  amount  of  bias  needed  for  photocurrent  readout. 

c.  The  stacks  of  each  band  should  be  limited  to  20  repetitions  of  the  basic 
barrier/well/barrier  cell  (due  to  the  complexity  to  fabricate  and  to  proc¬ 
ess  a  large  amount  of  different  crystal  layers). 

d.  The  stacks  should  be  grown  onto  the  substrate  in  the  following  se¬ 
quence:  NIR,  MWIR,  and  LWIR.  This  is  to  prevent  the  undesired  ab¬ 
sorption  mechanisms  from  sequentially  blocking  the  IR  bands  of  pass¬ 
ing  through  their  respective  detector  stacks.  Thus,  the  longer  bandgap 
stacks  should  be  placed  before  (bottom)  the  smaller  ones.  Also,  the 
MWIR  and  LWIR  should  be  the  second  (mid)  and  the  third  (top)  in  or¬ 
der  to  apply  the  respective  grating  pattern. 

e.  The  indium  compositions  in  InvGa/_<As  as  well  as  the  aluminum  compo¬ 
sitions  in  AlrGa/_vAs  should  be  restricted  to  no  more  than  3  different 
values  each  in  the  entire  crystal.  This  is  due  to  MBE  growing  limitations 
discussed  previously. 

f.  The  contact  layers  should  be  heavily  n-doped  (about  2. 1018  cm-3 )  in  or¬ 
der  to  improve  the  contact  uniformity  in  long  pixels. 

g.  The  cap  layer  (top)  and  the  contact  layer  between  the  MWIR  and  LWIR 
quantum  well  stacks  should  be  thick  enough  to  allow  the  application  of 
the  respective  grating  patterns. 

Following  the  considerations  given  above,  several  trials  were  executed  varying 
the  dimensions  and  the  compositions  of  each  band  separately.  The  final  stack  configura¬ 
tions  with  the  respective  energy  and  detection  wavelength  are  shown  in  Figures  17,  18, 
and  19.  Its  important  to  remark  that,  since  the  photodetectors  operate  cooled,  all  the  de¬ 
sign  computations  were  carried  using  77K  material  parameters. 
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Figure  17  Parameters  of  the  NIR  quantum  well  stack  (left)  and  its  respective  basic  cell 
potential  profile  with  the  confined  energy  levels  (right). 


The  NIR  quantum  wells  are  basically  asymmetric  step  wells  where  the  overlap  be¬ 
tween  the  initial  (valence  band)  and  final  (conduction  band)  envelope  functions,  the  ma¬ 
jor  contributor  for  the  strength  of  the  transition,  was  maximized,  bounded  by  the  limit  in 
indium  concentration  and  the  location  of  the  conduction  band  exited  state.  The  well  is  not 
doped  to  minimize  the  undesired  intersubband  transitions. 

The  MIR  quantum  wells  are  basically  symmetric  step  wells  where  the  intersub¬ 
band  oscillator  strength  was  maximized,  bounded  by  the  same  parameters  of  the  previous 
configuration.  In  this  case,  two  GaAs  layers  were  placed  on  the  sides  of  the  Ino.25Gao.75As 
layer  to  better  accommodate  the  AI0.40  Gao.6As/Ino.25Gao.75As  lattice  constant  mismatch 
and  to  reduce  the  amount  of  indium  in  the  structure.  All  the  three  layers  of  the  well  were 
heavily  doped  to  increase  the  absorption  coefficient. 
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Figure  18  Parameters  of  the  MWIR  quantum  well  stack  (left)  and  its  respective  basic 
cell  potential  profile  with  the  confined  energy  levels  (right). 
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Figure  19  Parameters  of  the  LWIR  quantum  well  stack  (left)  and  its  respective  basic  cell 

potential  profile  with  the  confined  energy  levels  (right). 


The  LWIR  quantum  wells  are  symmetric  and  shallow  square  wells  where  the  dop¬ 
ing  concentration  was  kept  lower  than  the  previous  configuration  to  reduce  the  leakage 
current  due  to  thermoionic  emission.  This  is  a  very  common  GaAs/AlGaAs  configuration 
[28]. 

The  entire  structure  is  composed  by  68  layers,  including  the  substrate.  Table  1 
compiles  the  parameters  needed  to  grow  the  crystal.  The  band  diagram  of  the  structure  is 
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shown  in  Figure  20.  The  allowed  transitions  for  each  quantum  well  configuration  are  also 
depicted  and  those  maximized  by  design  are  highlighted.  Notice  that,  in  the  figure,  the 
layer  numbers  correspond  to  those  in  the  Table  1. 


Table  1  GaAs  based  QWIP  sample  description. 


PROFILE 

THICKNESS 

A 

DOPANT 

CONC. 

±10% 

1 

n-GaAs 

7,000 

Si 

2xl018  cm'3 

2 

i-Alo.26Gao.74As 

300 

None 

3 

n-GaAs 

zYZU 

52 

Si 

0.5xl018  cm"3 

4 

i-Alo.26Gao.74As 

300 

None 

5 

n-GaAs 

5,000 

Si 

2xl018  cm'3 

6 

i-Alo.40Gao.60As 

300 

None 

7 

n-GaAs 

13 

Si 

2xl018  cm'3 

8 

n-Ino.15Gao.85  As 

zYZU 

14 

Si 

2x1 01S  cm'3 

9 

n-GaAs 

13 

Si 

2xl018  cm'3 

10 

i-Alo.40Gao.60As 

300 

None 

11 

n-GaAs 

5,000 

Si 

2xl0ls  cm'3 

12 

i-GaAs 

300 

None 

13 

i-Ino.10Gao.90  As 

X20 

43 

None 

14 

i-Ino.25Gao.75As 

40 

None 

15 

i-GaAs 

300 

None 

16 

n-GaAs 

7,000 

Si 

2xl018  cm'3 

17 
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Figure  20  Band  diagram  of  the  GaAs  based  QWIP  sample.  The  layer  numbers  corre¬ 
spond  to  those  described  in  Table  1. 
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Despite  the  restrictions  imposed  by  the  limitations  thoroughly  discussed,  the 
GaAs-based  sample  was  fabricated  by  means  of  contract  by  IQE  Inc.  and  succeeded  in 
covering  part  of  the  three  desired  spectral  regions  (the  results  are  presented  and  discussed 
in  the  next  chapter).  To  some  extent,  the  InP  system  presents  even  more  restrictions,  as 
described  in  the  following  section. 


2.  InP  System 

Semiconductor  structures  on  InP  substrates  are  more  difficult  to  grow,  fabricate, 
and  process.  Furthermore  the  band  offset  data,  one  of  the  most  important  parameter  in  de¬ 
signing  QWIPs,  is  very  controversial,  as  mentioned  before,  making  the  risk  in  such  at¬ 
tempts,  very  high.  Nonetheless,  the  bandgap  of  InP  is  smaller  than  GaAs,  making  this 
choice  very  attractive  for  the  NIR  region.  In  addition,  the  lattice  matched  combinations 
eliminate  the  strain  restrictions  widening  the  design  possibilities. 

Again  considering  the  discussed  possibilities,  restrictions,  and  the  lessons  learned 
with  the  GaAs  based  sample  results  (discussed  in  the  next  chapter)  the  following  design 
requirements  were  set  for  the  detectors  based  on  InP  system: 

a.  The  InP  based  QWIP  should  detect  only  two  IR  spectral  regions  in¬ 
stead  of  three  originally  set.  This  is  in  attempt  to  reduce  the  complex¬ 
ity  and  consequently  the  costs, 

b.  The  peak  absorption  of  the  two  different  bands  should  be:  above  1 .2 
pm  and  below  1.3  pm  for  NIR  (to  verify  the  ability  to  detect  higher 
wavelengths  in  this  band,  allowing  the  necessary  room  to  accommo¬ 
date  the  band  offset  uncertainties),  and  as  close  as  possible  to  the  mid¬ 
dle  of  the  band  for  MWIR  (~4.0  pm ). 

c.  For  the  MWIR  quantum  wells,  it  would  be  preferable  to  use  the  ternar¬ 
ies  In0.53Ga0.47As  (well)  and  Al0.45ln0.52As  (barrier)  due  to  the  consis¬ 
tency  in  their  properties  data  [20]. 

d.  Bound-to-continuum  transitions  should  be  preferred  to  broaden  detec¬ 
tion  spectrum  of  the  MWIR  detection. 


64 


e.  The  stacks  of  each  band  should  have  no  more  than  40  and  no  less  than 
30  repetitions  of  the  basic  barrier/well/barrier  cell.  This  is  to  increase 
the  absorption  coefficient  peak,  naturally  lower  than  the  previous  de¬ 
sign,  caused  by  the  bound-to-continuum  transitions. 

f.  The  stacks  should  be  grown  onto  the  substrate  in  the  following  se¬ 
quence,  NIR  and  then  MWIR.  Beyond  the  reasons  mentioned  in  the 
previous  design,  the  MWIR  has  to  be  on  the  top  of  the  device  to  re¬ 
ceive  the  grating  pattern. 

g.  The  contact  layers  doping  concentration  should  be  reduced 

(1.1018  cm-3)  in  contrast  to  the  GaAs-based  QWIP,  to  reduce  the  prob¬ 
ability  of  electrons  diffusing  to  lower  potentials,  populating  the  non- 
doped  wells.  This  was  experienced  in  the  GaAs-based  sample.  Also 
the  bottom  layer  should  be  p-doped  in  order  to  create  a  p-i-n  structure 
and  deplete  of  electrons  the  NIR  wells  reducing  the  possibility  of  un¬ 
desired  intersubband  transitions  even  more. 

h.  A  stop-etch  layer  should  be  included  right  below  the  cap  layer  to  im¬ 
prove  the  precision  in  processing  the  grating  pattern.  This  was  not  con¬ 
sidered  in  the  previous  design  counting  on  the  higher  controllability  on 
etching  GaAs-based  semiconductors. 

Attempting  to  consider  all  the  defined  requirements,  several  trials  (using  the  com¬ 
putational  tool)  were  executed  varying  the  dimensions  and  the  compositions  of  each  band 
independently.  The  final  stack  configurations  are  shown  in  Figures  21  and  22.  The  tem¬ 
perature  considerations  are  the  same  as  the  previous  design. 

Similar  to  the  previous  design,  the  NIR  quantum  wells  are  basically  asymmetric 
step  wells.  In  this  case,  there  are  no  strain  limitations;  the  asymmetry  is  clearly  higher 
than  the  GaAs  sample,  favoring  the  desired  interband  transitions.  Again,  the  well  is  not 
doped  to  minimize  the  undesired  intersubband  transitions. 
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Figure  21  Parameters  of  the  NIR  quantum  well  stack  (left)  and  its  respective  basic  cell 
potential  profile  with  the  confined  energy  levels  (right). 
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Figure  22  Parameters  of  the  MWIR  quantum  well  stack  (left)  and  its  respective  basic 
cell  potential  profile  with  the  confined  energy  levels  (side). 

The  MWIR  quantum  wells  are  highly  doped,  symmetric,  and  deep  wells.  In  this 
case  there  is  only  one  confined  level  forcing  the  bound-to-continuum  transitions.  Notice 
that  the  ternaries  In0.53Ga0.47As  (well)  and  Af  ^In^As  (barriers)  provide  an  absorption 
peak  very  close  to  the  middle  of  the  band.  Similar  configurations  were  reported  success¬ 
fully  in  [32], 
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The  entire  structure  is  composed  by  77  layers,  including  the  substrate.  Table  2 
compiles  the  parameters  needed  to  grow  the  crystal.  The  band  diagram  of  the  structure  is 
shown  in  Figure  23.  The  allowed  transitions  for  each  quantum  well  configuration  are  also 
depicted,  and  those  maximized  by  design  are  highlighted.  Notice  that  in  the  figure,  the 
layer  numbers  correspond  to  those  in  the  Table  2. 


Table  2  InP  based  QWU 

P  sample  description. 

PROFILE 

THICKNESS 

A 

DOPANT 

CONC. 

±10% 

1 

n-InP 

3,800 

Si 

lxl 018  cm"3 

2 

n-Alo.4gIno.56As 

300 

Si 

lxlO18  cm"3 

3 

i-Al0  48In0.56As 

X40 

300 

None 

4 

n-Ino.53Gao.47As 

35 

Si 

2x1 018  cm'3 

5 

i-Alo.48lno.56As 

300 

None 

6 

n-InP 

7,000 

Si 

lxlO18  cm'3 

7 

i-  (Ino.53Ga0.47As)0.5o(InP)  0.50 

300 

None 

8 

i-  (Ino.53Ga0.47As)0.7o(InP)  0.30 

X30 

85 

None 

9 

i-ino.53Gao.47As 

30 

None 

10 

i-  (Ino.53Ga0.47As)0.5o(InP)  0.50 

300 

None 

11 

p-InP 

7,000 

Be 

lxlO18  cm"3 

12 

S.I.  InP  Substrate 

-1st  QW stack-  tp  -  4.2  rm- 
- 40  x - J 


-2Ta  QW  stack  -?p~  1.3  rrm- 
- 30  x - J 


Conduction  Band 


4.2^im 


0.85|nm 


26.2|jm 


1.33|nm 


1.18|nm 


1.27(im 


10 


Valence  Banc 

11  12 


Figure  23  Band  diagram  of  the  InP  based  QWIP  sample.  The  layer  numbers  correspond 

to  those  described  in  Table  2. 
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Although  the  InP  sample  is  more  attractive  to  meet  the  requirements,  the  imple¬ 
mentation  of  such  crystalline  structure  is  very  expensive  and  risky,  due  to  the  reasons 
previously  discussed.  In  an  attempt  to  make  the  photodetector  feasible,  it  was  restricted  to 
a  two-color  detector  and  data  as  reliable  as  possible  were  used  from  [20],  [32],  and  [45], 
In  addition,  the  feedback  from  the  GaAs  sample  absorption  measurements  induced  some 
changes  in  the  design  approach,  making  the  InP  a  rich  tool  of  learning.  The  sample  was 
sent  to  fabrication  and  its  characterization,  testing  and  processing  are  matter  for  future 
work. 


This  chapter  covered  the  aspects  of  the  QWIPs’  design.  First,  groups  of  semicon¬ 
ductors  were  analyzed  and  strained  GaAs-based  and  matched  InP-based  systems  proved 
more  suitable  to  simultaneously  detect  the  NIR,  MWIR  and  LWIR.  Heterostructures 
formed  by  stacks  of  uncoupled  wells  were  chosen  and  the  sequence  of  detection  bands 
was  defined.  Second,  considering  the  specificities  of  the  selected  configurations,  compu¬ 
tational  tools  were  developed  and  the  structure  self-consistent  Schrodinger-Poisson  solu¬ 
tions  were  obtained  by  the  shooting  method.  Finally,  two  QWIP  samples  were  designed. 
The  first,  for  three-color  detection,  is  based  on  GaAs/AlxGa;.xAs/  InvGa/-.vAs  on  GaAs 
substrate  and  the  second,  for  two-color  detection,  is  based  on  ( Inft 5 jG aft 4 7 A s)v( In P) j_y/ 

Alo. usIno^As  on  InP  substrate.  The  GaAs-based  sample  was  fabricated  via  contract  by 
IQE  Inc.  and  its  characterization  is  addressed  in  the  next  chapter  along  with  the  design  of 
the  photodetector  devices. 
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IV.  EXPERIMENTAL  RESULTS  AND  FABRICATION 


The  GaAs  QWIP  sample  was  grown  by  way  of  contract  by  IQE  Inc.  using  MBE 
on  a  semiinsulating  GaAs  wafer  with  3  inches  diameter,  following  the  specifications  dis¬ 
cussed  in  the  last  chapter.  The  room  temperature  absorption  measurements  were  carried 
out  and  the  results  are  discussed  in  this  chapter  along  with  the  optical  coupling  and  the 
mask  design  for  detector  processing. 


A.  ROOM  TEMPERATURE  ABSORPTION  MEASUREMENTS 

There  are  many  ways  to  characterize  QWIP  samples.  One  way  is  to  measure  the 
transitions  that  can  occur  inside  the  crystal.  As  shown  in  Figure  20,  the  GaAs  sample  ex¬ 
hibits  both  interband  and  intersubband  possible  transitions.  Cathodeluminescence  and 
photo  luminescence  [12]  can  be  used  to  measure  interband  transitions  between  ground 
states  and  Fourier  transform  infrared  spectroscopy  (FTIR)  [11]  can  be  used  to  measure 
interband  transitions.  Very  comprehensible  descriptions  of  the  FTIR  working  principle 
are  given  in  [11]  and  [53],  and  can  be  summarized  as  follows. 

The  infrared  absorption  measurements  of  the  MIWR  and  LWIR  bands  were  car¬ 
ried  out  in  SRL  using  a  Nicolet-750  FTIR.  Unlike  normal  optical  spectroscopy,  which 
uses  diffraction  grating  monochromators  in  the  visible  and  ultraviolet,  FTIR  uses  an 
Michaelson  interferometer.  A  broadband  radiation  source  is  directed  towards  a  beam 
splitter  where  the  light  is  divided  into  two  equal  beams.  One  beam  is  directed  to  a  fixed 
mirror  and  the  other  to  a  moving  mirror.  The  beams  are  reflected  back  to  the  beam  split¬ 
ter  which  combines  them  and  directs  them  towards  the  sample.  A  detector,  placed  behind 
the  sample,  measures  the  combined  beam  intensity  after  passing  through  the  sample.  [53] 

In  this  configuration,  the  detector  sees  the  original  signal  superimposed  by  a  time- 
delayed  copy  of  itself.  This  represents  the  autocorrelation  function  of  the  radiation  field, 
commonly  called  interferogram.  The  spectrum  is  obtained  by  taking  the  Fourier  trans¬ 
form  of  the  interferogram.  The  center  of  the  interferogram  corresponds  to  the  white  light 
position  (WLP),  where  the  optical  path  is  the  same  in  both  branches  of  the  Michaelson  in¬ 
terferometer.  The  entire  source  radiation  band  is  reflected  towards  the  sample.  The  phase 
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difference  between  partial  waves  outside  the  WLP  is  wavelength  dependent.  The  spectral 
resolution  of  the  FTIR  is  determined  by  the  displacement  length  of  the  moving  mirror. 
[53] 

The  FTIR  spectrometer  uses  a  computer  to  control  the  movable  mirror  and  record 
the  detector  data.  The  mirror  is  moved  through  its  entire  range  at  a  speed  controlled  by 
the  user  and  the  detected  intensity  is  recorded  as  a  function  of  data  points.  The  data  points 
are  converted  to  a  time  scale  using  the  velocity  of  the  mirror  and  the  sampling  frequency. 
This  procedure  is  repeated  several  times  and  the  data  are  averaged  to  improve  the  signal- 
to-noise  ratio.  The  controlling  software  performs  the  Fourier  transform  of  the  resulting 
interferogram.  [11] 

In  order  to  obtain  only  the  quantum  wells’  absorption,  the  undesired  effects  must 
be  removed  from  the  measured  interferogram.  That  includes  the  atmospheric  and  bulk 
material  absorptions.  The  easiest  way  to  do  that  is  to  use  a  polarizer  in  front  of  the  detec¬ 
tor.  By  setting  the  polarization  direction  parallel  to  the  quantum  well  layers,  no  absorp¬ 
tion  is  expected  from  the  quantum  wells.  Then,  all  that  is  measured  is  considered  back¬ 
ground  and  subtracted  from  the  actual  measurement. 

The  samples  were  prepared  in  the  form  of  a  4.79  mm  long,  13.6  mm  wide  and 
0.66  mm  thick  waveguide  by  polishing  45°  angle  facets  to  enhance  the  absorption,  fol¬ 
lowing  the  same  procedures  described  in  [1 1],  The  fabricated  sample  waveguide  is  shown 
in  Figure  24. 

The  prepared  sample  is  placed  on  a  gold-plated  holder  that  channels  the  infrared 
beam  from  spectrometer  into  the  waveguide.  The  light  enters  the  device  normal  to  the  45° 
face,  reflects  through  the  waveguide  and  transmits  out  the  opposing  face  (as  depicted  in 
Figure  24)  towards  the  polarizer  and  detector. 
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Figure  24  GaAs  sample  waveguide  prepared  for  the  absorption  measurements.  (Sche¬ 
matic  diagram  on  the  top  and  the  actual  sample  on  a  holder  on  the  bottom. 


The  absorptance  of  the  sample  was  measured  at  room  temperature  at  polarizations 
from  zero  (TM)  to  90°  (TE)  at  steps  of  10°.  The  results  are  presented  in  Figure  25. 


Wavelength  (urn) 

Figure  25  MWIR  and  LWIR  absorptance  of  the  GaAs  based  sample  measured  at  90°  in¬ 
cident  polarizations.  The  absorption  peaks  are  highlighted. 
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The  absorption  peaks  for  MWIR  and  LWIR  are  5.275  pm  and  8.745  pm,  respec¬ 
tively.  Notice  that  the  presence  of  an  additional  high  LWIR  absorption  with  peak  around 
13.8  pm.  The  presence  of  this  absorption  is  discussed  later  in  this  section.  The  depend¬ 
ence  of  the  absorption  peaks  strength  with  the  polarization  angle  is  shown  in  Figure  26. 
The  absorption  peak  strength  decreases  with  the  square  of  cosine  of  the  polarization  an¬ 
gle,  confirming  that  the  absorptions  are  due  to  the  quantum  wells  (discussion  in  Section 
II.A.6). 


Figure  26  MWIR  (a)  and  LWIR  (b  and  c)  absorption  peak  strengths  as  a  function  of  the 

polarization  angle. 


At  first  view,  some  interesting  facts  can  be  identified.  The  MWIR  absorption  at 
5.3  pm  was  expected  to  be  stronger  than  the  LWIR  due  to  its  higher  doping  concentra¬ 
tion.  The  LWIR  peak  at  13.8  pm  was  not  expected  since  the  quantum  wells  responsible 
for  that  absorption  are  not  doped  (see  Table  1  and  Figure  20).  A  possible  explanation  for 
these  observations  is  summarized  in  Figure  27. 

The  figure  shows  the  configuration  of  the  sample  layers  with  their  doping  concen¬ 
trations  and  the  band  diagram.  The  transitions  of  interest  are  represented  by  black  arrows. 
The  graphical  insertion  represents  the  measured  absorption  in  the  sample.  The  straight 
red  arrows  indicate  which  well  stack  is  responsible  for  which  absorption  region.  Also  the 
unexpected  LWIR  absorption  peak  is  highlighted  along  with  its  intersubband  transition. 
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Figure  27  Schematic  diagram  of  the  effects  of  the  electrons  drift  and  tunneling  to  the 
NIR  and  LWIR,  respectively,  due  to  the  high  doping  concentration  of  the  con¬ 
tact  layers. 

Due  to  the  difference  in  potential  between  layers,  the  electrons  from  the  highly 
doped  surrounding  contact  layers  can  diffuse  into  the  NIR  quantum  wells  (indicated  in 
Figure  27  by  the  “curved”  double  red  arrows),  populating  these  wells  with  electrons,  al¬ 
lowing  the  intersubband  transition  around  13.9  pm .  Also,  the  electrons  can  tunnel 
through  the  AlGaAs  barriers  to  the  LWIR  quantum  wells  (indicated  in  Figure  27  by  the 
“curved”  single  red  arrows),  increasing  the  carriers’  populations  in  those  layers.  With  all 
quantum  wells  having  nearly  the  same  electron  concentration,  the  strength  of  the  transi¬ 
tions  is  directly  proportional  to  the  thickness  of  the  layer  and  the  transition  wavelength. 
This  explains  the  enhancement  of  absorption  peak  strength  with  increasing  wavelength. 

Those  effects  are  undesired  since  they  will  contribute  to  increase  in  the  dark  cur¬ 
rent  due  to  the  high  probability  of  thermoionic  emission  from  shallow  wells.  One  ap¬ 
proach  to  avoid  that  in  future  designs  is  to  reduce  the  contact  layers  doping  concentration. 
In  addition  a  p-doping  in  the  bottom  layer  would  make  a  p-i-n  structure  which  depletes 
electrons  in  the  NIR  quantum  wells,  eliminating  the  probability  of  intersubband  transition 
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there.  This  approach  was  used  in  the  InP-based  QWIP  design,  as  discussed  in  Section 
III.E.2. 


The  absorption  coefficient,  a ,  can  be  obtained,  for  each  band,  given  that  the  ab- 
sorptance,  A,  is  [11] 

A  =  -  log„(e-“1'),  (4.1) 

where  L  is  the  path  lengh  through  the  quantum  wells,  represented  by 

Lwg  -  nNL cos(45°) .  (4.2) 

Here,  n  is  the  number  of  quantum  wells  in  each  pass,  N  is  the  number  of  passes  and  L  is 
the  well  period.  The  absorption  coefficient  of  the  sample  was  computed  and  compared 
with  the  theoretical  calculations  as  shown  in  Figure  28. 


Figure  28  Comparison  between  the  estimated  and  the  measured  MWIR  and  two  LWIR 

absorption  coefficients. 


The  blue  curve  is  the  estimated  absorption  using  Equation  (2.66)  and 
r  =  20  meV .  The  red  curve  is  the  estimated  absorption  considering  the  effects  of  the  dif¬ 
fusion  and  tunneling  discussed  above  as  well  as  using  the  measured  broadening  parame¬ 
ters.  The  uncertainties  in  the  estimation  of  the  absorption  coefficient  are  several;  how¬ 
ever,  it  is  possible  to  have  an  idea  of  the  general  behavior.  Notice  that  in  Figure  28  the 

red  curve  is  shifted  from  the  black  by  almost  the  same  amount  in  all  bands. 
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On  the  other  hand,  the  estimated  peak  positions  are  within  0.3  pm  of  the  meas¬ 
ured  values,  which  indicates  the  model  gives  good  description  of  the  absorption.  A  better 
visualization  is  possible,  comparing  the  normalized  curves,  as  shown  in  Figure  29. 


Figure  29  Comparison  between  the  estimated  and  the  measured  MWIR  and  LWIR  nor¬ 
malized  absorption. 


The  FTIR  measurement  did  not  cover  the  NIR  region  due  to  the  limitations  on  the 
equipment.  However,  the  detection  in  this  band  can  be  characterized  after  the  devices’ 
fabrication  through  the  responsivity  measurement. 

The  absorptance  measurement  was  important  to  confirm  the  MWIR  and  LWIR 
detection  bands.  Also  the  differences  between  the  estimated  and  the  measured  absorption 
peaks  can  be  used  to  refine  the  low  temperature  material  parameters  and  design  the  opti¬ 
mized  light  coupling  gratings.  The  recomputed  low  temperature  absorption  peaks  are  5.0 
pm  and  8.5  pm .  This  is  important  because  the  grating  acts  as  a  narrow  bandpass  filter 
and,  if  not  well  tuned  with  the  absorption  peak  wavelength,  it  can  degrade  the 
photodetector  performance.  The  optical  coupling  design  is  addressed  in  the  following 
section. 
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B.  OPTICAL  COUPLING 

As  mentioned  in  Section  II.A.6  and  confirmed  by  the  results  showed  in  the  previ¬ 
ous  section,  the  selection  rules  state  that  intersubband  absorption  cannot  occur  in  sym¬ 
metric  structures,  for  TE  mode  incident  radiation.  In  practical  devices,  detection  of  nor¬ 
mal  incident  IR  is  desired.  In  this  design  it  is  required.  The  most  efficient  way,  that  has 
been  reported  [29]  [30]  [31]  to  overcome  this  selection  rule,  consists  in  etching  diffrac¬ 
tion  gratings  on  the  top  of  the  device.  Several  configurations  have  being  reported  and  the 
most  common  are:  random  grating  [29],  when  pattern  uniformity  is  the  main  goal  (impor¬ 
tant  for  arrays);  corrugated  grating  [30],  when  multi-band  detection  is  the  priority;  and  bi¬ 
directional  ordered  cross  grating  [50],  when  one  band  detection  is  desired.  The  later  con¬ 
figuration  is  simpler  than  the  others,  in  terms  of  fabrication,  and  it  has  shown  potentiali¬ 
ties  for  multi-band  detection  [55],  Therefore,  bi-dimensional  ordered  grating  was  chosen 
for  our  GaAs-based  QWIPs,  to  be  implemented  in  two  different  patterns,  one  maximized 
for  LWIR  and  other  maximized  for  MWIR. 

A  thorough  theoretical  analysis  is  given  in  [55],  It  is  concluded  there  that  the  cou¬ 
pling  efficiency  is  very  sensitive  to  the  grating  period  and  the  dimensions  of  the  cavities. 
Several  fabrication  limitations,  such  as  smallest  possible  dimension  and  precision  in  etch¬ 
ing,  compel  us  to  employ  a  non-optimum  grating  with  a  simpler  architecture  [29]. 

Although  the  two  grating  parameters  are  dependent  on  each  other  [55],  they  are 
treated  separately  here  without  disturbing  significantly  the  final  result  [29],  First  the  grat¬ 
ing  period  is  analyzed,  followed  by  the  area  and  depth  of  the  groves.  The  parameters  are 
optimized  for  the  detection  bands  and  adapted  to  meet  the  fabrication  limitations. 

A  schematic  diagram  of  the  diffraction  grating  in  a  quantum  well  device  is  shown 
in  Figure  30.  There,  the  substrate,  contact  layers,  a  quantum  well  stack  and  a  metal  layer 
are  sketched  exactly  as  they  are  placed  in  actual  devices.  The  dimensions  are  off-scale  to 
better  show  the  beam  paths. 
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Figure  30  Schematic  diagram  of  the  diffraction  grating 


In  the  figure  S  represents  the  groove  width,  h  represents  the  groove  depth,  and 
A  represents  the  period  of  the  grating.  The  backdiffracted  beam  angle  of  the  zth  order  is 
given  by  the  classical  grating  equation  [56] 

sin  6®  +  sin  6t  =  n  — —  .  (4.3) 

nr  A 


Here  6i  is  the  incident  beam  angle  with  respect  to  the  normal,  n  is  the  diffraction  order, 

and  nr  is  the  refraction  index  of  the  material.  Consequently,  for  normal  incidence  the  dif¬ 
fracted  angles  exhibit  the  following  dependence  on  the  grating  period: 


0D  =  arcsin 


V 


A 

nA 


(4.4) 


According  to  Equation  (4.4),  the  maximum  first-order  broadside  diffraction  oc¬ 
curs  when  the  grating  period  has  the  same  size  of  the  incident  wavelength  inside  the  ma¬ 
terial.  Hence,  the  optimum  value  of  A  is 


A 


op 


A 

nr  ' 


(4.5) 


Figure  3 1  shows  the  first-  and  second-order  diffraction  angle  dependence  on  the  grating 
period,  normalized  by  the  wavelength  inside  the  material. 
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Figure  3 1  Variation  of  the  first  and  second  order  diffracted  angles  with  the  normalized 

grating  period. 

The  high  sensitivity  of  the  diffraction  angles  with  the  variation  of  the  grating  pe¬ 
riod  is  clearly  shown  in  Figure  31,  mainly  for  periods  smaller  than  1.5  A.  This  becomes 
important  when  the  desired  precision  in  fabrication  cannot  be  obtained.  Using 
nGraAs  -  3.21  in  Equation  (4.5),  the  optimum  grating  period  for  MWIR  (peak  at  5.0  pm ) 
and  LWIR  (peak  at  8.5  pm)  were  determined  as  1.5  and  2.6  pm ,  respectively.  The  next 
step  is  to  determine  the  area  and  depth  of  the  groove. 

The  refraction  grating  can  be  viewed  as  two  surfaces  where  one  is  retarded  rela¬ 
tive  to  the  other  by  a  distance  h,  the  groove  depth.  If  one  surface  area  is  u  times  smaller 
than  the  other,  the  area  duty  cycle  turns  out  to  be  1/(1  +  u ) .  Using  the  same  classical 
analysis  used  to  compute  the  interference  pattern  of  a  Fabry-Perot  etalon  [56],  the  optical 
difference  path,  Ag ,  between  two  beams  reflected  by  the  two  surfaces  can  be  expressed 

as 

Ag  =  nrh(  1  +  cos  6)  (4.6) 

where  9  is  the  angle  between  the  diffracted  beam  and  the  normal  (see  Figure  30).  Con¬ 
sequently,  the  electric  field,  proportional  to  the  sum  of  the  two  beams  can  be  written  as 
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(4.7) 


c  ik(x+ A  )  ikx 

tg  oc  e  s  +ue 

where  k  is  the  wave  vector  and  x  is  the  arbitrary  distance  in  the  direction  of  k.  The  inten¬ 
sity  is  proportional  to  FgFg  *  and  can  be  represented  by 

Ig  oc  (1  +  u2)  +  2ucos(kAg) .  (4.8) 

An  optimum  grating  is  the  one  that  minimizes  the  back  reflection  ( I  — »  0,  when  6^  — >  0 ) 
and  maximizes  the  broadside  reflection  {6  — »  90°) . 

On  one  hand,  when  0  -  0 ,  the  optical  path  difference  reduces  to 

Ag  =  2 nrh  ,  (4.9) 

and  Equation  (4.8)  can  be  rewritten  as 

^  +  M  -  =  -  cos(2 knrh)  when  0  =  0° .  (4.10) 

2  u 

Since  u  represents  a  positive  quantity  and  the  left  side  of  Equation  (4.10)  must  be  in  the 
interval  (-1,1) ,  the  only  possible  value  for  u  is  one  (uop  =  1 ),  indicating  that  the  projected 

area  of  the  cavity  openings  should  be  one  half  of  the  total  projecting  area.  Therefore,  it  is 
possible  to  determine  the  optimum  value  of  h.  This  result  is  expressed  as 

hop=^  +  n^~  n  =  0,  1,2,3,...  .  (4.11) 

4  nr  2  nr 

Figure  32  shows  the  back-reflected  intensity  due  to  grating  effects.  The  parameters  u  and 
h  are  normalized  by  their  deduced  optimum  values. 
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Figure  32  Intensity  back-reflection  due  to  diffraction  grating. 

It  is  rather  interesting  to  observe  the  behavior  of  the  backreflected  intensity  when  h  is 
varied.  Figure  33  shows  this  dependence  for  the  MWIR  and  LWIR  detection  peaks  along 
with  the  GaAs  QWIP  layer  structure. 


Figure  33  Dependence  of  the  grating  backreflected  intensity  on  the  cavity  depth.  The 
GaAs  QWIP  layers  were  placed  right  below  the  graph  to  compare  the  depth 

relative  to  the  actual  device. 
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In  Figure  33,  the  optimum  points  for  LWIR  detection  are  points  3  and  7  and  for 
MIWR  detection  are  1,4  and  6,  indicated  by  arrows.  Points  2  and  5  are  the  best  choices  if 
a  single  grating  would  be  applied  for  both  bands.  Due  to  fabrication  limitations,  the  grat¬ 
ing  must  be  etched  from  the  top  of  the  wafer,  while  the  sample  is  flat  (i.e.,  before  etch¬ 
ing).  Consequently,  for  LWIR  the  h  value  corresponding  to  point  1  should  be  selected. 

For  MWIR,  point  1  is  intuitively  the  best  choice;  however,  to  be  able  to  read  the  photo¬ 
current  for  each  band  separately  as  indicated  in  Figure  7,  the  grating  cavity  must  reach 
the  second  (from  the  top)  contact  layer  in  order  to  make  contact  and  allow  an  independent 
readout.  As  a  result,  the  only  choice  that  maximizes  LWIR  is  an  h  corresponding  to  point 
6.  This  situation  is  depicted  schematically  in  Figure  34. 


Figure  34  Schematic  diagram  of  the  grating  configuration  to  allow  three  separate  photo¬ 
current  readouts.  The  vertical  is  exaggerated  for  visual  purposes. 

An  important  factor  that  must  be  considered  is  the  etching  tolerance.  Since  in  both 
points,  3  and  6  are  very  close  to  the  respective  quantum  well  stacks.  Therefore  a  small  re¬ 
duction  in  the  grating  deepness  should  be  considered  in  order  to  assure  that  the  cavity 
will  not  reach  the  quantum  wells,  even  though  the  grating  efficiency  degrades.  The  opti¬ 
mum  and  the  selected  values  of  h  are  in  Table  3  in  the  end  of  this  section. 


81 


To  complete  the  grating  specifications,  the  horizontal  dimension  of  the  cavity,  as 
well  as  the  distance  between  them  must  be  computed.  Figure  35  shows  a  schematic  dia¬ 
gram  of  grating  pattern.  There,  Al  and  A2  represent  the  area  to  receive  the  grating  pattern 

and  the  groove  area,  respectively,  and  x  is  the  distance  between  grating  squares,  called 
the  critical  dimension  (CD). 


□ 


□ 


Figure  35  Schematic  diagram  of  a  bi-dimensional  ordered  grating  pattern. 


Considering  that  u  =  1 ,  the  relationship  between  the  areas  Ax  and  A2 ,  depicted 
in  Figure  35  is 

A,  =2(NM)A2.  (4.12) 

Using  Equation(4.12)  and  after  some  algebra,  the  expressions  for  the  dependence  of  the 
remaining  parameters  with  relation  to  the  grating  period  is  given  by 


and 


(4.13) 


x  =  0.29A . 


(4.14) 
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Having  the  grating  period  previously  calculated,  the  remaining  parameters  were  com¬ 
puted  and  the  results  were  summarized  in  Table  3. 

A  restriction  to  the  grating  design  is  the  minimum  dimension  that  can  be  achieved 
during  fabrication  stage.  Based  on  the  limitations  discussed  in  the  next  section,  the  mini¬ 
mum  feature  allowed  in  this  design  was  set  as  1  pm .  The  grating  horizontal  parameters 
were  recalculated  for  x  =  1  pm  and  the  results  are  presented  in  Table  3.  This  imposes  dif¬ 
ferent  losses  in  efficiency  for  each  detection  band  as  can  be  predicted  using  Figure  33. 


Table  3  Grating  parameters  of  the  GaAs  sample.  (All  units  are  pm) 


Optimum 

Feasible  (tolerance  =  ±0.1  pm ) 

A 

V 

X 

h 

A 

V 

X 

h 

MWIR 

1.5 

1.1 

0.45 

1.91 

3.4 

2.4 

1.0 

1.8 

LWIR 

2.6 

1.8 

0.75 

0.65 

3.4 

2.4 

1.0 

0.6 

The  grating  design  was  an  effort  to  allow  the  normal  incidence  detection  for  the 
QWIP  devices.  It  was  not  an  isolated  task,  since  it  is  very  dependent  on  the  fabrication 
processes  and  the  configuration  of  the  devices.  Actually  an  iterative  work  was  necessary 
to  adjust  the  grating  to  the  devices  and  the  devices  to  the  grating  fabrication  possibilities. 
Although  the  parameters  that  maximize  the  detection  on  each  band  were  degraded  from 
their  optimum  values  to  make  the  devices  feasible,  normal  incidence  detection  still  can  be 
achieved  for  both  bands. 

C.  DEVICES  FABRICATION 

The  ultimate  goal,  after  extensive  work  was  to  fabricate  the  photodetectors  in  or¬ 
der  to  execute  the  performance  analysis  for  each  band  separately.  The  device  fabrication 
sequence,  procedures,  rules,  and  limitations  must  be  considered  at  the  very  first  steps  of 
heterostructure  design.  Initially,  the  knowledge  of  the  processing  possibilities  and  con¬ 
straints  can  determine  the  thickness  of  the  contact  layers,  the  addition  of  stop-etch  layers 
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in  the  structure,  and  even  the  sequence  of  the  quantum  well  stacks.  Lastly,  the  devices’ 
shapes  and  configurations  are  strongly  dependent  on  the  processing  issues. 


Our  devices  will  be  fabricated  by  the  National  Research  Council  in  Canada  and 
their  characterization  is  left  for  future  work  due  to  lack  of  time;  however  the  photodetec¬ 
tor  configurations  were  defined  and  the  mask  layouts  were  drawn,  based  on  the  following 
requirements: 

a.  The  mask  plates  should  have  the  dimensions  of  be  4x4  inches  (limited 
by  the  mask  aligner). 

b.  The  useful  area  of  the  masks  cannot  be  grater  than  44x26  mm  (limited 
by  the  aligner) 

c.  The  smallest  feature  permitted  in  the  mask  is  1.0  pm  (limited  by  the 
mask  fabrication  costs  and  plate  dimensions). 

d.  The  grating  patterns  should  be  applied  on  the  top  of  the  wafer,  while  the 
surface  is  still  flat  enough.  Consequently  the  grating  etching  should  be 
the  first  step. 

e.  The  devices  should  be  grouped  in  not  more  than  4x4  mm  (limited  by  the 
package). 

f.  The  contact  areas  should  be  not  less  than  50x50  pm .  This  is  a  limitation 
imposed  by  the  wire  bonding  equipment. 


In  addition,  considering  the  necessity  of  three  separate  readouts  and  2  grating  pat¬ 
terns,  the  general  processing  sequence  was  defined  [57],  simulated  using  MEMS  Pro  v5.0 
licensed  for  NPS.  The  results  are  shown  in  Figure  36. 
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Wafer  after  MBE  (68  layers) 
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Coating  with  photoresist 
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Transfer  the  second  step  mesas  pattern 

3 

Transferring  the  first  grating  pattern 

16 

Etching  the  second  step  mesas 

4 

Etching  the  first  grating 

17 

Removing  the  remaining  photoresist 
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Removing  the  remaining  photoresist 
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Coating  with  photoresist 
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Transfer  the  third  step  mesas  pattern 
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Transferring  the  second  grating  pattern 
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Etching  the  third  step  mesas 
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Etching  the  second  grating 
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Removing  the  remaining  photoresist 
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Removing  the  remaining  photoresist 
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Coating  with  photoresist 
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Coating  with  photoresist 
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Transfer  the  top  mesas  pattern 
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Etching  the  top  mesas 
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Metal  liftoff 

Figure  36  Photodetector  devices  fabrication  sequence. 


Finally,  to  be  able  to  test  several  properties  of  the  QWIPs  such  as  three  color  de¬ 
tection,  normal  incidence  detection,  effects  of  contacts  in  long  pixels,  detection  uniform¬ 
ity,  effect  of  the  bias  in  the  detection  band,  grating  pattern  effect  on  the  other  band,  and 
polarization  and  angle  of  incidence  effects  during  the  characterization  phase,  different 
photodetectors  configuration  were  designed.  Figure  37  shows  the  basic  detectors. 
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For  practical  reasons,  the  fabrication  was  divided  in  two  phases.  In  the  first  phase, 
devices  without  gratings  will  be  fabricated  and  in  the  second  phase,  after  the  preliminary 
tests,  the  devices  with  a  grating  pattern  will  be  fabricated.  However,  to  fit  in  the  research 
budget,  only  two  masks  accommodating  all  layers  necessary  for  both  phases  will  be  fab¬ 
ricated.  The  first  mask  is  clear  field  with  the  features  in  chrome,  for  the  mesas  etching. 
The  second  mask  is  dark  field  with  the  features  clear,  for  the  grating  patterns  and  metal 
layers. 

The  characterization  of  the  devices  will  be  carried  out  by  SRL  and  the  results  will 
be  published  in  future  work. 


This  chapter  presented  the  characterization  of  the  GaAs-based  QWIP  sample, 


showing  absorption  on  the  desired  bands  and  some  effects  of  the  doping  and  band  struc- 
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ture  of  the  sample.  The  results  were  important  as  feedback  for  the  next  designs.  Having 
the  structures  tested  and  working,  the  optical  coupling  was  designed  for  the  optimum 
normal  incidence  detection.  Limitations  on  the  fabrication  forced  the  final  patterns  to  be 
away  from  the  optimum  points,  but  not  compromising  the  overall  proof-of-concept.  Fi¬ 
nally  the  photodetector  devices  were  designed  to  permit  the  measurement  of  the  figures 
of  merit  and  several  other  effects.  Despite  the  limitations,  the  design  methodology  de¬ 
fined  during  this  research  and  the  experimental  results  show  a  great  potentiality  of  this 
approach,  addressed  in  the  conclusion  chapter. 
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V.  CONCLUSIONS 


The  realm  of  QWIP  design  is  very  extensive  [28].  The  advances  in  material  sci¬ 
ences  and  semiconductor  growth  technology  have  allowed  the  use  of  bandgap  engineer¬ 
ing  to  cover,  ideally,  all  IR  bands  [5].  Specifically  in  the  area  of  QWIPs  a  great  effort  was 
made  in  the  past  20  years  to  make  those  devices  commercially  available  for  several  appli¬ 
cations.  High  sensitivity,  high  selectivity,  and  multispectral  detection  capability  are  some 
of  the  properties  of  quantum  well  detectors  that  make  those  devices  very  attractive  for 
military  applications.  The  possibility  of  combine  NIR,  MWIR  and  LWIR  in  a  single  de¬ 
tector,  motivated  the  work  described  in  this  documented  in  this  text.  The  successful  ac¬ 
complishments  along  with  the  lessons  learned  and  the  future  possibilities  are  discussed  in 
this  chapter. 


A.  FINAL  CONSIDERATIONS 

The  purpose  of  this  thesis  was  to  study,  to  model,  and  to  design  a  QWIP  capable 
to  detect  three  different  IR  bands  within  the  wavelengths  intervals  of  0. 9-1.4  pm ,  3. 8-5.0 
pm  and  8.0-12.0  pm .  This  objective  was  accomplished  through  the  following  steps. 

In  the  first  chapter,  the  background  theory  was  organized  in  order  do  physically 
explain  the  quantum  effects  in  semiconductor  heterostructures  that  make  possible  the  IR 
detection.  The  envelope  function  approximation  was  used  to  mathematically  model  the 
electron  behavior  under  the  influence  of  the  band-offset  and  charge  distribution  poten¬ 
tials.  The  one-dimensional  Schrodinger’s  equation  was  used  to  obtain  the  confined  en¬ 
ergy  levels  and  their  respective  wavefunctions.  The  transitions  between  different  states 
due  to  the  incident  IR  radiation  were  identified  and  studied.  The  absorption  coefficient 
was  modeled  for  interband  and  intersubband  transitions,  and  the  strength  and  selection 
rules  were  identified.  The  effects  of  the  electric  field  were  also  included  in  the  mathe¬ 
matical  models.  The  mechanisms  of  extraction  of  the  electrons  from  the  wells  and  leads 
to  a  photocurrent  were  also  studied  and  modeled.  After  that,  the  main  figures  of  merit, 
necessary  to  set  requirements  and  to  analyze  the  QWIPs  performance  were  discussed. 
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The  second  chapter  brings  the  practical  design  implications  and  procedures.  The 
semiconductor  material  properties  along  with  their  availability  were  considered  and  two 
main  material  groups  were  selected  to  be  used.  GaAs/AlGaAs/InGaAs  system  showed 
easy  in  growth  and  fabrication  while  InP/AlGaAs/AlInAs  system  was  capable  of  meeting 
the  requirements  for  detecting  the  desired  NIR  band.  A  configuration  of  three  uncoupled 
quantum  wells  stacks,  one  stack  for  each  band,  placed  between  highly  doped  contact  lay¬ 
ers,  showed  to  be  more  suitable  to  allow  separate  readouts  of  the  signals. 

Considering  the  specificities  of  the  configuration  and  materials,  a  computational 
tool  was  developed  to  solve  self-consistently  the  Schodinger-Poisson  equations  for  the 
structure.  The  shooting  method  was  implemented  due  to  its  ability  to  handle  any  potential 
profile,  necessary  feature  for  QWIP  designs.  Finally,  a  GaAs  based  sample  for  three  color 
detection,  and  a  InP  based  sample  for  two  color  detection,  were  designed,  taking  into  ac¬ 
count  the  limitations  on  each  configuration. 

The  GaAs  sample  was  fabricated  using  MBE  and  the  absorption  measurements 
were  presented  in  the  fourth  chapter.  An  FTIR  spectrometer  was  used  to  measure  the 
room  temperature  absorptance  of  the  MWIR  and  LWIR  bands,  and  the  sample  showed 
absorption  in  good  agreement  with  the  design  predictions.  Also  the  heavily  doped  contact 
layers  introduced  a  second  LWIR  absorption  band.  The  measured  peaks,  5.3  pm ,  8.7  pm 
and  13.8  pm  varied  less  than  0.3  pm  from  the  predicted  values  mostly  due  to  the  uncer¬ 
tainties  in  the  material  parameters. 

The  optical  coupling  structure  was  then  designed  to  allow  detection  of  the  nor¬ 
mally  incident  radiation.  A  bi-dimensional  ordered  grating  pattern  was  selected  and  opti¬ 
mized  separately  for  both  MWIR  and  LWIR  desired  peaks.  Constraints  in  mask  fabrica¬ 
tion  and  wafer  processing  forced  the  degradation  of  the  grating  performance,  to  make  it 
feasible.  Finally  the  photodetector  device  configurations  were  designed,  to  permit  to  per¬ 
formance  analysis.  The  fabrication  and  characterization  of  the  prototypes  are  matter  for 
future  work. 

Despite  the  financial,  temporal  and  technological  limitations,  normally  present  in 
any  research,  the  objective  of  this  thesis  was  successfully  accomplished.  However,  fur¬ 
ther  work  needs  to  be  done  to  finalize  the  devices’  fabrication  and  testing  and  to  improve 
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the  mathematical  models  and  the  numerical  algorithms  to  be  able  develop  optimal  detec¬ 
tors.  This  is  summarized  as  follows. 

B.  FUTURE  WORK 

The  natural  follow-on  task  of  this  research  is  the  GaAs-based  devices’  fabrication. 
The  designed  configurations  allow  the  analysis  of  the  NIR  absorption,  through  the  photo¬ 
current  spectroscopy.  In  addition,  figures  of  merit  can  be  analyzed  along  with  several  ef¬ 
fects  such  as:  MWIR  and  LWIR  grating  efficiency,  MWIR  grating  influence  in  LWIR  de¬ 
tection,  LWIR  grating  influence  in  MWIR  detection,  uniformity  of  grating  patterns,  inci¬ 
dent  radiation  polarization  effects,  direction  of  incidence  effects,  etc.  All  the  results  can 
be  used  to  optimize  the  design  methodology  and  models  to  improve  the  design  accuracy. 

Another  important  aspect  is  the  characterization  of  the  InP-based  sample.  Since 
that  design  incorporates  modifications  induced  by  the  feedback  of  the  GaAs  sample 
measurements,  this  would  be  the  natural  way  to  test  and  validate  the  new  approach. 

In  addition,  the  models’  weakness  and  points  of  improvement  identified  in  the 
text  can  be  summarized  as  following: 

a.  The  numerical  computations  were  done  considering  the  valence  non¬ 
degenerate  with  respect  to  the  heavy  and  light  holes’  bands.  The  sim¬ 
plest  way  to  include  the  degeneracy  effects  is  to  solve  numerically  the 
4x4  k»p  Hamiltonian  described  in  Equation  (2.16)  to  obtain  the  valence 
band  confined  levels. 

b.  The  heterostructure  electrostatic  potential  was  computed  without  ac¬ 
counting  for  the  strain  effects.  This  can  be  included  in  the  model  follow¬ 
ing  the  approach  discussed  in  [53], 

c.  The  absorption  coefficient  for  interband  transitions  were  computed  only 
for  the  TE  mode.  The  TM  absorption  can  be  included  using  a  more  rig¬ 
orous  approach  suggested  in  [24], 
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d.  The  Stark  shift,  included  in  the  absorption  coefficient  algorithms  was 
not  properly  verified  against  experimental  results.  This  can  be  done 
along  with  the  GaAs-based  devices’  responsivity  measurements. 

e.  The  photocurrent  as  well  as  the  figures  of  merit  estimation  should  be  in¬ 
cluded  in  the  computational  tools.  The  mathematical  models  can  be  re¬ 
fined  after  the  GaAs-based  devices  characterization. 


Finally,  the  work  done  through  the  models,  computational  tools  and  design  ap¬ 
proach  discussed  in  this  thesis,  corroborated  by  the  experimental  results,  can  be  itera¬ 
tively  improved  by  fabrication  and  characterization  of  new  devices. 
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